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1.  INTRODUCTION 


— 

J-  [This  report  describes  an  improved  linear  analytical  model  and 
digital  program  developed  for  the  calculation  of  axisymmetric  ^launch 
vehicle  steady-stage  response  to  applied  axisymmetric  sinusoidal  loads.  / 

- — -  L  s 

The  detailed  computer  programming  manual  for  the  digital  program  is 
contained  in  Volume  II  of  this  report. 


In  the  evaluation  of  launch  vehicle  behavior,  it  is  necessary  to  study 

the  response  of  the  entire  vehicle  to  a  wide  variety  of  dynamic  loadings  to 

insure  the  structural  integrity  and  stability  of  the  system.  Much  effort 

has  already  gone  into  the  development  of  techniques  to  calculate  the  vehicle 

response  to  lateral  and  longitudinal  loadings  using  distributed  and  lumped 

spring -mass  models  and  techniques  for  theoretical  and  empirical  modeling 

12 

of  the  vehicle  behavior.  9  However,  experimental  data  indicate  that 
these  procedures  are  unsatisfactory  in  several  respects.  For  example, 
accurate  representation  of  important  structural  shell  characteristics  and 
realistic  coupling  of  the  fluids  with  the  detailed  structural  behavior  of 
tank  walls  and  bulkheads  are  omitted. 


The  approach  described  herein  overcomes  the  above  noted  defi¬ 
ciencies.  A  finite  element  technique  is  utilized  to  construct  the  total 
launch  vehicle  stiffness  matrix  £kJ  and  mass  matrix  JmJ  by  subdividing 
the  prototype  structure  into  a  set  of  (1)  axisymmetric  shell  components, 
(2)  fluid  components,  and  (3)  spring-mass  components.  In  this  way,  it 
is  possible  to  represent  as  separate  shell  units  the  fairing,  interstage 
structure,  bulkheads,  tank  walls  and  engine  thrust  structure,  and  to  con 
veniently  provide  for  the  inertial  and  stiffness  characteristics  of  equip 
ment,  engines  and  vehicle  supporting  structure 

r**-  ‘7* 

4fj  dtec 


The  stiffness  and  mass  matrices“Tor  the  complete  launch  vehicle 
are  obtained  by  superposition  of  the  stiffness  and  mass  matrices  of  the 
individual  shell,  fluid  and  spring -mass  components  which  are  computed 
using  a  Rayleigh-Ritz  approach.  Fluid  motions  are  assumed  to  be  con¬ 
sistent  with  the  shell  component  distortions.  The  superposition  technique 
assures  displacement  compatibility  and  force  equilibrium  at  the  joints 
between  components.  After  the  complete  system  stiffness  and  mass 
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matrices  have  been  formulated,  displacement  boundary  conditions  are 
introduced  by  removing  appropriate  rows  and  columns  corresponding  to 
points  on  the  vehicle  and  its  supports  which  are  rigidly  restrained  from 
motion. 

The  coupled  system  natural  frequencies  and  mode  shapes  are 
obtained  from  the  eigenvalue  equation  constructed  with  the  total  stiff¬ 
ness  and  mass  matrices 


[k]  {a}  -  p2  [m]  (a)  =  0  (1.1) 

in  which  p  is  the  circular  frequency  of  the  system  and  is  the  modal 

vector  whose  components  are  the  longitudinal,  radial  and  rotational  dis¬ 
placements  at  discrete  points  on  the  vehicle.  The  steady-state  response 
due  to  simple  harmonic  loads  is  determined  using  a  standard  modal  re¬ 
sponse  procedure  which  expresses  the  total  displacement,  velocity, 
acceleration  and  force  responses  as  the  linear  superposition  of  the  in¬ 
dividual  modal  responses  based  on  an  assumed  modal  damping. 

The  procedure  will  handle  shell  components  with  a  wide  range  of 
geometries.  It  includes  shell  effects  in  the  tank  and  bulkhead  structure, 
but  avoids  the  need  for  including  detailed  local  deformation,  such  as  at 
shell  discontinuities,  which  are  unimportant  in  determining  the  total 
dynamic  behavior  of  the  vehicle.  The  approach  has  the  capability  of  rep¬ 
resenting  the  tank  or  stage  of  most  interest  in  great  detail  and  those  of 
least  interest  with  minimum  detail,  as  desired,  thereby  minimizing  the 
computation  time  required  and  remaining  within  the  maximum  limitations 
of  standard  eigenvalue  routines.  The  formulation  of  the  problem  is  sub¬ 
divided  into  well-defined  portions,  leading  to  efficient  coding  and  easy 
modification  for  later  incorporation  of  asymmetric  shell  behavior  and 
even  more  detailed  treatment  of  the  fluid  behavior. 

The  analytical  procedure  discussed  herein  is  summarized  in  Figure  1. 
The  launch  vehicle  analytical  model  is  discussed  in  Section  2.  The  equations 
for  the  shell  and  fluid  component  stiffness  and  mass  matrices  are  developed 
in  Sections  3  and  4,  respectively.  The  coordinate  representation  which 
forms  the  framework  for  the  vehicle  model,  and  the  construction  of  stiff¬ 
ness  and  mass  matrices  for  the  complete  launch  vehicle,  are  discussed 


in  Sections  5  and  6,  respectively.  The  method  for  computing  the  dynamic 
response  is  discussed  in  Section  7.  The  computer  program  arrangement 
is  described  in  Section  8  and  detailed  input  data  requirements  for  the  com¬ 
puter  program  are  itemized  in  Section  9.  A  complete  list  of  symbols  is 
provided  in  Section  10. 


Figure  1.  Summary  of  Technical  Approach 


1.4 


2.  ANALYTICAL  MODEL 


As  illustrated  in  Figure  2,  the  vehicle  structure  is  subdivided  into 
a  consistent  set  of  shell  components,  a,  fluid  components,  b,  and  multi - 
coordinate  spring-mass  components,  c.  The  total  vehicle  which  may  be 
represented  is  limited  to  one  with  not  more  than  six  (6)  fluid  components. 
The  total  number  of  shell  components  which  may  be  represented  shall  not 
exceed  forty  (40).  The  characteristics  of  the  spring -mass  components 
representing  such  equipment  as  engines  and  mass-elastic  supports  are 
provided  directly  by  low  order  (<  10)  stiffness  and  mass  matrices.  The 
total  number  of  spring -mass  components  may  not  exceed  thirty  (30). 

The  vehicle  behavior  is  described  in  terms  of  motions  of  discrete  points 
on  the  vehicle  located  at  intersections  of  shell  components,  at  lumped 
masses,  and  at  intermediate  points  on  the  shell  elements.  The  number 
of  nonfixed  degrees -of -freedom  by  which  the  behavior  of  the  system  is 
described  may  not  exceed  eighty  (80). 

The  specific  shell  components  to  be  used  are  conical  frustums 
(which  include  cylindrical  shells  as  a  special  case)  and  ellipsoidal  bulk¬ 
heads  (which  include  hemispherical  shells  as  a  special  case).  Within 
the  domain  of  thin  shell  theory,  the  shell  components  may  have  ortho¬ 
tropic  properties  and  a  linear  thickness  variation  in  the  meridional 
direction.  Local  thickening  of  a  shell  at  a  bulkhead  or  wall  joint  may  be 
handled  by  using  an  equivalent  local  hoop  stiffener  which  is  provided  as 
input  in  the  form  of  an  additional  spring -mass  component.  Initial  static 
stresses  based  on  membrane  theory  are  accounted  for  in  determining  the 
stiffness  matrix  for  the  shell  components. 

The  most  general  fluid  component  may  be  in  contact  with  an  ellip¬ 
soidal  upper  bulkhead,  a  conical  tank  wall,  and  a  conical  or  ellipsoidal 
lower  bulkhead.  The  bulkhead  shell  elements  may  be  convex  down  or  up 
with  fluid  at  any  desired  depth  on  either  side,  both  sides  or  neither  side. 
The  tank  configurations  which  are  considered  are  illustrated  in  Figure  3. 
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B.  SYSTEM  COMPONENTS 


Figure  2.  Vehicle  Components 
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Figure  3.  Tank  Configurations 


3.  SHELL  COMPONENT  STIFFNESS  AND  MASS  MATRICES 


The  displacements  of  an  individual  shell  component  (Figure  4)  are 
approximated  in  the  Rayleigh-Ritz  manner  by  a  finite  series  of  functions 
having  the  form 


U  _ 

u(g)  =  £  *  u  <  11 

k=i  -  K 

(3.  1) 

V  __ 

V(g)  =  £  P,v  U)  ,  v<  11 

i=l 


in  which  £  is  a  dimensionless  variable  and  0  <  1.  For  conical  frustums 

Figure  5,  £  =  s I JL  sin  For  convex  upward  ellipsoidal  shells,  Figure  6, 

£  =  < p/<pQ>  and  for  convex  downward  ellipsoidal  shells,  £  =  (tt  -  0)/(tt  -  0^). 

The  assumed  mode  shapes  u^(^)  and  v^(^)  consist  of  polynomial 
terms  sufficient  to  represent  all  modes  of  shell  distortion,  including 
longitudinal  stretching,  radial  dilatation  and  rigid  body  displacements. 

The  specific  shape  of  the  assumed  modes  is  determined  within  the 
limitations  of  a  tenth  order  polynomial,  as  follows: 
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Uk(^  =  L  akn£n 


n=0 


10 


(3.2) 


v<(5>=  £0bin«n 

n=0 


in  which  [aj  =  [a](1Jx11)  and  [b^J  =  [b](VxU)  define  the  poly- 
nomial  functions  associated  with  the  local  coordinates  (a^.)  and  (|3^), 
respectively.  TaJ  and  are  furnished  as  input  to  provide  maximum 

flexibility  in  the  selection  of  the  assumed  coordinate  functions.  It  should 
be  understood  that  the  a^,  and  the  (3^  are  the  unknown  generalized  coordi¬ 
nates  associated  with  the  shell  component,  whereas  the  a^n  and  b ^  are 
definite  arbitrarily  specified  coefficients  which  determine  the  functions 


3.  1  SHELL  STIFFNESS  MATRIX 


The  shell  stiffness  matrix  is  constructed  from  the  potential  energy 
function  V  derived  in  Appendix  A.  For  the  generalized  displacement 
(Equation  3.  1),  the  stiffness  matrix  is  defined  by 


((U+V)x(U+V)) 
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which  is  consistent  with  the  Rayleigh-Ritz  procedure  and  in  which,  for 
example 
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and 


k  =  1  .  .  .  U 
i  =1  ...  V 

The  strains  (e  ,  gQ),  curvatures  (K  ,  K  )  and  meridional  rotation  p  are 
(p  U  <p  u 

defined  in  terms  of  the  displacements  in  Appendix  A  by  Equations  (A.  13) 
through  (A.  17).  The  quantities  C^,  C^,  C^,  C^,  ^34  an^  ^44  are 
the  orthotropic  stress -strain  coefficients  which  are  functions  of  the 
dimensionless  parameter  For  the  present  analytical  model,  these 
coefficients  are  approximated  by  the  following  polynomial  equations: 

'cnW 

{5(D)  =  /  C12(g)  \  ={C(0)}  +  i  ({c<l>}  -  (c<0)})  (3.5) 
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The 

are  provided  as  input 

N^°  is  the  initial  meridional  stress  caused  by  tank  pressures  and 
longitudinal  vehicle  accelerations  (prestress  and  static  longitudinal 
acceleration  stress).  These  initial  stresses  are  derived  in  Appendix  B 
using  membrane  theory,  which  is  a  reasonable  first  order  approximation, 
except  in  very  localized  areas  where  the  bending  stresses  predominate. 
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This  simplifies  the  calculation  of  the  initial  stresses  which  are  readily 
obtained  from  the  membrane  equations  of  equilibrium.  Expressions 
for  the  initial  stresses  for  conical  shell  components  are  given  in 
Equations  (  B.  1)  through  (B.  6),  for  the  upright  ellipsoidal  bulkheads  in 
Equations  (  B.  7)  through  (B.  11),  and  for  the  inverted  ellipsoidal  bulkheads 
in  Equations  (B.  12)  through  (B.  15).  The  computation  of  these  stresses 
is  included  as  an  integral  part  of  the  digital  program. 

After  substitution  of  the  displacements  (3.  1)  into  (3.  3),  the  shell 
component  stiffness  matrix  becomes 

[kJ  =  2* Jr  (kl  •  {u}-  (u}T  +  K2  (u)  •  (u}T  +  K3  [{u}*  (ii)T  +  (u)  •  {u}T] 

4-  K4  j^v}  •  (v}T  4-  {V}  *  (v}T  j  +  K5  ^v^  •  ^v^  4-  K6 

4-  K7  ^v}  •  {v}T  4-  K8  J^v}  •  (v}T  4-  ^v}  •  (v}T  j 

4-  K9  [{v}  •  {u}T  4-  ^u^  •  (v}T]  +  KlO^v^  •  (u}T  +  ^u^  • 

+  K  a  [{V>  •  (i.)T  +  (i)  •  (5}T]  +  K 1 2  [{<-}  •  (u}T  +  (u)  •  {v}T] 

+  K13[{5}-  {u}T  +  {u}  •  {»}T]')r1d0 

L  '  (3.6) 

in  which 
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for  a  cone. 


The  constants  K1  -  K13  are  defined  as  follows: 


2  if  *1  \ 

K1  =  Cu  sin  <p  +  C  ^  — 2  ^sin  0  +  —  cos  (pj 


*34  rr. 


sin  0  +  —  cos  0  +  C 

ri  7 


COS  0  ,  TVT  o  2 

44 - 2  +  N0  cos 

r  r 


K2  =  C  cosL  0 


K3  =  -C33  ^sin  0  +  ^  cos  <pj  +  C. 


K4  =  C 


12  r 


K5  = C22  ~J 

r 


Z  i  /  ri  \ 

K6  =  ^  cos  0  +  C33  — ^  (cos  0  -  —  sin  0  j 

r,  V  rl  / 


34  rr, 


r,  \  •  Z  A 

cos  0  -  —  sin  0  I  +  C44 - 

1  /  r 


+  N  °  sin^  0 
0 


K7  =  C  ~  ~  sin  0 


k»  -  c33  ^  (~*  -  £  *»  *)  +  c34  sial 


K9  =  -C10 

12  r 


3.6 


K10  -  -C  j  j  sin  0  cos  0  -  j  |cos  0  -  sin  0^  ^sin  0  +  cos  0^ 

•  ^  ,  1  A  ,  ~  cos  0  f  rl  .  A 

sm  0  +  —  cos  0j  +  C34  — (cos  0  -  —  sin  0j 


,  -  sin  0  cos  0  ,  - T  o  . 

+  c-44 - E— 2 -  +  N  sin  0  cos  0 

r  ^ 


K11  =  _C33  ^T“  ^sin  0  +  ~  cos  0^  sin  0  +  C34  sin-^  ^°S — & 


K12  =  C„  fcos  0  -  li  sin  ^  +  r  sin  ft  H°s2  fr 

33  ^  r  r  y  r.  34  r 

K  13  =  C33  sin  0  cos  0  (3.8) 

3.  2  SHELL  MASS  MATRIX 

The  shell  mass  matrix  associated  with  the  generalized  coordinates 
ak’  ^ i  *s  derived  by  operating  on  the  expression  for  the  kinetic  energy  T. 
F or  each  shell  component,  the  mass  matrix  is  defined  by 


3.7 


where  y  is  the  shell  mass  density. 

The  detailed  equations  for  the  calculation  of  the  stiffness  and  mass 

matrices  of  the  shell  components  are  presented  in  Appendix  D.  The 

equations  are  written  in  matrix  notation  to  accommodate  the  Gaussian 
5 

weighted  matrix  integration  scheme. 


X 


Figure  4.  Displacements  of  Shell  of  Revolution 


x 


LONGITUDINAL 

UPPER  EDGE  OF  VEHICLE  AXIS 

SHELL  COMPONENT 


(TRANSVERSE  CROSS-SECTION  THROUGH  LONGITUDINAL  AXIS) 

Figure  5.  Conical  Shell  Component 
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LONGITUDINAL 


(TRANSVERSE  CROSS-SECTION  THROUGH  LONGITUDINAL  AXIS) 
Figure  6.  Ellipsoidal  Bulkhead  Component 
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4.  FLUID  COMPONENT  MASS  MATRIX 


In  addition  to  the  mass  and  stiffness  matrices  for  the  shell  compo¬ 
nents,  the  inertial  effects  due  to  the  presence  of  liquid  propellants  in  the 
vehicle  fuel  tanks  must  be  considered.  The  linear  analytical  model  does 
not  include,  however,  the  effective  fluid  stiffness  caused  by  changes  in 
the  fluid  head  during  shell  distortions  as  this  is  a  higher  order  nonlinear 
effect. 

The  general  fluid  component  b  is  enclosed  by  three  shell  components, 
consisting  of  conical  and  ellipsoidal  shells  of  revolution.  The  specific  tank 
configurations  which  are  included  in  the  present  model  have  been  illustrated 
in  Figure  3.  For  a  typical  fluid  component,  as  shown  in  Figure  7,  the  tank 
is  divided  into  three  shells  in  which  the  upper  bulkhead  is  referred  to  as 
shell  al,  the  tank  wall  as  shell  a 2,  and  the  lower  bulkhead  as  shell  a3. 

The  fluid  motions  are  a  function  of  the  generalized  displacements 
for  shell  components,  al,  a2  and  a3.  For  the  or  v^(£)  °f  eac^  shell 

element  defined  by  Equation  (3. 2),  there  is  an  associated  fluid  motion 
u  (x)  and  v  (x,  r),  where  l.miW  and  W  =  (U  +  V)  ,  +  (U  +  V)  ->  + 

(U  +  V)  g-  ^m(x)  fluid  displacement  parallel  to  the  x-axis  (longitudi¬ 

nal  axis  of  the  vehicle)  associated  with  the  tank  shell  generalized  displace¬ 
ment  m.  vm(x,  r)  is  the  fluid  displacement  parallel  to  the  r-axis  (radial 
axis  of  the  vehicle)  associated  with  the  tank  shell  generalized  displace¬ 
ment  m. 

The  general  form  of  the  fluid  mass  matrix  consistent  with  the  fluid 

3 

displacement  may  be  expressed  as 


is  the  fluid  mass  density  and  matrix  °f  order  W. 


The  fluid  motion  um(x)  is  assumed  independent  of  location  r  and  is 
obtained  by  treating  the  fluid  as  incompressible  and  inviscid.  From  these 
assumptions,  ^m(x)  is  equal  to  the  change  in  volume  below  a  given  loca- 
tion  x  divided  by  the  corresponding  tank  cross  sectional  area.  Thus  (see 
Figure  7), 

bjdx  (4.3) 

where 


{^x0  =  -  -4  f  _n  r(cot  *  {u<^}b+  {v(^} 


Fluid  sloshing  motions  which  disturb  the  planar  character  of  the 
assumed  longitudinal  motion  are  beyond  the  scope  of  this  treatment,  but 
may  be  superimposed  as  generalized  sloshing  modes  independent  of  the 
shell  distortions.  Consistent  with  the  assumed  longitudinal  fluid  motion 

u,  (x)  and  the  axisymmetric  nature  of  the  linear  model,  the  radial  fluid 

K  .  A  .  Av  ,  .  ,  h 

motion  v^.(x,  r)  varies  linearly  with  space  coordinate  r.  At  a  particular 

longitudinal  location  x,  the  radial  fluid  motion  is  a  function  of  the  radial 

motion  of  the  adjacent  tank  shell  boundary  and  of  the  longitudinal  fluid 

motion.  Thus 

{v(x,  r)}  =  J  ^cot  <p  {u(£)}b  +  {v(£))b  -  cot  ?>{u(x)}j  (4.5) 


Upon  substituting  Equations  (4.3)  and  (4.5)  into  (4.1)  and  integrating 
with  respect  to  r,  one  obtains 


(4.6) 


where 

=  -  J_-  r  (cot  *  C<F>}b  +  {^>b  )dx 

and 

{rv(x,  r iy  =  r  cot  0  {  u(£)}b  +  "  r  cot  0  {u(x)}  (4.  7) 

The  detailed  expressions  developed  for  evaluating  Equation  (4.6)  for 
various  cases  are  summarized  in  Appendix  C. 


SHELL  1 

(UPPER  BULKHEAD) 


Figure  7.  Definition  of  Flu: 
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5.  LAUNCH  VEHICLE  COORDINATE  SYSTEMS 


5.  1  LOCAL  COORDINATE  DISTORTIONS 

The  matrices  for  the  shell  and  fluid  components  are  derived  in  Sec¬ 
tions  3  and  4,  respectively,  using  a  system  of  generalized  coordinate  dis¬ 
placements  (local  coordinates),  as  given  by  Equations  (3.  1)  and  (3.2)  to 
describe  the  shell  distortion.  These  equations  may  be  written  in  matrix 
notation  in  the  following  form: 

T 


where 


u<e)=  [akJ  [A]{e}=[£j[A]{aJ 

(5.1) 

we)  =  |pJ  H (0  =  UJ  [b ]T{Pi> 

(5.2) 

|akJ  =  [5X  •  •  •  a-J 

|/ij =  L^i  *  ■  •  M 

(5.3) 

[sj  =  [i(e)(D2.  •  •  (e)10J 

U  and  V  are  constants  provided  as  input  which  define  the  number  of  local 
coordinates  selected  for  representing  shell  distortions  in  the  longitudinal 
and  radial  directions,  respectively,  for  a  particular  component. 

5.2  SYSTEM  COORDINATE  POINT  DISPLACEMENTS 

In  order  to  work  with  reference  to  a  space  frame,  however,  it  is 
necessary  to  transform  the  generalized  coordinates  to  space  coordinates 
designated  as  system  coordinates.  The  system  coordinates  represent  dis¬ 
placement  and  rotations  at  specific  points  on  each  shell  component,  con¬ 
nections  of  spring-mass  components,  and  applications  of  force  inputs  at 
arbitrary  stations  along  the  vehicle. 

For  each  shell  component,  system  coordinates  are  provided  to  repre¬ 
sent  displacements  in  the  longitudinal  and  radial  directions  at  equally  spaced 
intervals  along  the  shell  meridian,  and  tangential  rotations  at  the  edges  of 
the  shell  (see  Figures  8,  9  and  10).  Longitudinal  displacements  (u^i0 


at  locations  =  (U  -  i)/ (U  -  1),  where  i  =  1,  2,  .  .  . ,  U  are  expressed  fusing 
Equation  (5.  1 ) j  as 


where 


(u(ii,V*u=  mm  (O 


w- 


q  q  *  •  ■  (q ) 


(5.4) 


(5.5) 


Ju  ^U*  ’  ■ 


Longitudinal  displacements  Av(|  )^  at  locations  £.  =  (V  -  j)/(V  -  1), 

J  J 

where  j  =  1,  2,  .  ..,  V,  are  similarly  expressed  (using  Equation  (5.2)  as 


where 


(v(eAv*if  MO]  OJ 


q  q  •  •  •  (q) 


(5.6) 


[v]  = 


(5.7) 


l'v  ev  •  •  •  (ev)  J 

The  scalar  quantities  U  and  V  are  constants  provided  as  input  which 
define  the  number  of  system  coordinate  point  displacements  to  be  provided 
in  the  longitudinal  and  radial  directions,  respectively,  for  a  given  shell 
component.  These  are  related  to  U  and  V,  as  discussed  in  Section  5.3. 
The  shell  rotations,  pi  and  p£,  are  evaluated  at  £  =  1  and  0,  respectively, 
and  by  Equation  (A.  17)  can  be  shown  to  have  the  following  form: 


{P!}=  NOfOO*  KlOf  if} 


(5.8) 


where 


9  10 


0  0 


(2x11) 


and 


9  10 


0  0 


(2x11) 


D 


1 


cos  0  sin  0 
ro  ro 


r 


for  conical  frustums 


(5.9) 


COS  0 

—  for  convex  upward  ellipsoidal  bulkheads,  and 
o 


cos  0 
r  o 

for  convex  downward  ellipsoidal  bulkheads 


sin  0 

=  - j - -  for  conical  frustums, 

sin  0 
r  o 

—  for  convex  upward  elipsoidal  bulkheads,  and 
o 


sin  0 

7 — r - /  •  Q - r  for  convex  downward  ellipsoidal  bulkheads. 

(ri)  %  - n) 

K  1  o 

where  ^r^j  is  r^  evaluated  at  0  =  0q 


The  total  vector  of  system  coordinates  displacements  for  shell  com¬ 
ponent  "a"  is  defined  from  Equations  (5.4),  5.  6)  and  5.  8)  as 
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(5.10) 


As  indicated  in  Figures  8,  9  and  10,  each  displacement  and  rotation 
coordinate  is  identified  by  a  number  C  where  1  2S  C  2  and  N^.  is  defined 
as  the  total  number  of  system  coordinates  used  in  the  vehicle  model.  The 


identification  numbers  C  may  be  arbitrarily  specified  on  the  vehicle  inde¬ 
pendent  of  location  and,  hence,  must  be  identified  with  the  vector  ^ 
components  in  the  input  data.  For  this  purpose,  a  vector  is  provided 

for  each  shell  component  a  with  the  identification  numbers  illustrated  in 


Figures  8,  9  and  10  arranged  consistent  with  Equations  (5.4),  (5.6)  and 
(5.  10).  Thus,  in  general, 


For  ellipsoidal  bulkhead  components,  vector  elements  C^y  and 
C  are  omitted. 

p2 

5.  3  LOCAL  TO  SYSTEM  COORDINATE  TRANSFORMATION 
FOR  SHELL  COMPONENT 

The  local  coordinate  displacements  for  shell  component  Mau  are  re¬ 
lated  to  the  system  coordinate  displacements  by  the  transformation  £tJ 
as  follows: 

Consolidating  the  arrangement  of  local  coordinate  displacements,  as 
defined  in  Equations  (5.  1),  (5.  2)  and  (5.  3),  let 


5.4 


(5.12) 


(a}-_  ■ 

&y(U+V)  x  1 

Consistent  with  the  system  coordinate  displacement  vector  <“a)  given  by 
Equations  (5.4),  (5.6)  and  (5.8),  the  transformation  matrix  JtJ  which 
relates  local  to  system  coordinates,  is  defined  by  the  equation 


{“a}  =  W'Oa} 


(5.13) 


where,  in  general, 


[Ta]' 


[u] 


w 


_o_ 

M 

W 


w 


W 


(5.14) 


Thus 


{“a)  *  [Ta]{°a> 


(5.15) 


It  is  apparent  that  computation  of  the  matrix  T  requires  inversion 

r  >  i  L  aJ 

of  the  matrix  I  T  I  ,  which  must  therefore  be  nonsingular,  square,  and 

of  order  (U  +  V).  This  may  require  modification  of  the  general  Equation 

(5.14),  as  discussed  below.  Special  attention  must  be  given  the  scalar 

quantities  U,  V,  U,  V  and  the  matrices  £bJ,  ^D^Jto  SdLtisfy  the 

above  conditions  consistent  with  the  shell  component  boundary  conditions. 

For  ellipsoidal  bulkheads,  the  ^A^and^Bj  matrices  are  selected  so 
that  the  displacement  v(0)  and  the  rotation  p2  are  equal  to  zero  to  satisfy 
the  conditions  imposed  by  axial  symmetry.  This  implies  that  for  ellip¬ 
soidal  bulkheads  U  =  U,  V  =  V,  which  requires  that  the  last  row  of  j" vj, 

|dJ  and  Equati-on  (5.  14)  be  removed.  For  a  conical  bulkhead,  (bJ 

is  selected  so  that  v(0)  is  zero.  This  case  requires  that  the  last  row  of 
jV]  be  removed,  and  U  +  1  =  U,  V  =  V.  For  cylindrical  shells,  U  =  U, 

V  =  V  -  2.  For  other  cases,  U  +  V  +  2  =  U  +  Visa  necessary  condition. 
The  following  example  matrices  fulfill  the  above  requirements. 


For  nonbulkhead  components: 


10000000000 


J(Uxll) 


0  10  0 


0  0  10. 


(5.16) 


0  0  0  1 


(Illustrated  for  U  =  4) 


HlV: 


(Vxl  1 ) 


10  0. 


0  10  0 


0  0  1  0  0. 


0  0  0  1  0  0. 


0  0  0  0  1  0. 


0  0  0  0  0  1  0 


(5.17) 


(Illustrated  for  V  =  6) 


For  ellipsoidal  bulkhead  components: 


1  0  0 


l(Uxll) 


0  0  10 


0  0  0  1 


(5. 18) 


0  0  0  0  1 


(Illustrated  for  U  =  4) 


0  10  0 


[bU 


(Vxl 1 ) 


0  0  10 


0  0  0  1  0 


(5.19) 


0  0  0  0  1 


(Illustrated  for  V  =  4) 
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For  the  simple  unit  diagonal  M  and  [bJ  matrices  illustrated,  the 
matrix  |tJ  may  be  poorly  conditioned  and  difficult  to  invert  accurately 
if  U  and/or  V  are  equal  to  or  greater  than  six  (6).  For  these  cases  it  is 
recommended  that  Shifted  Chebyshev  Polynomial  Coefficients^  (see 
Table  1)  be  used  for  the  £aJ  and^B^  matrices  to  improve  the  accuracy  of 
the  matrix  |tJ  calculation.  Using  Table  1,  Equations  (5.16),  (5.17), 

(5.  18)  and  (5.  19)  become,  respectively, 

For  nonbulkhead  components: 


^(Uxll) 


fBJ(Vxll)  “ 


o 

o 

o 

i— H 

I _ 

0  0 

0  0 

-12  0  0 

• 

• 

1-880 

•  * 

•  • 

_-l  18  -48  32 

• 

•  • 

(Illustrated 

for  U  = 

4) 

0 

0 

• 

•  • 

2 

0  0 

• 

•  • 

-8 

8  0 

0 

• 

18 

-48  32 

0 

0 

-32 

160  -256 

128 

0 

50 

-400  1120 

-1280 

512  0 

(Illustrated  for  V  =  6) 
For  ellipsoidal  bulkhead  components: 


(5.20) 


.  .  0 


(5.21) 


M(Uxll)  = 


10  0  .  . 

1  0  8  0. 

-1  0  -48  32  0 

1  0  160  -256  128 
(Illustrated  for  U  =  4) 


(5.22) 


0  2 


0 


0 


0 


[bU 


(Vxll) 


0-8  8  0 
0  18  -48  32  0 

0  -32  160  -256  128 
(Illustrated  for  V  =  4) 


(5.23) 


0 


5.4  LOCAL  TO  SYSTEM  COORDINATE  TRANSFORMATION 
FOR  FLUID  COMPONENT 

As  illustrated  in  Figure  7,  and  discussed  in  Section  4,  the  distor¬ 
tion  of  each  fluid  component  nbn  is  a  function  of  the  distortion  of  the  three 
enclosing  shell  components  al,  a2  and  a3.  The  local  to  system  coordinate 
transformation  matrix  jT^J  for  the  fluid  component  is  thus  obtained  as  a 
combination  of  the  shell  component  transformations  for  a  =  al,  a2  and 

a3,  as  follows: 

Consolidating  the  arrangement  of  local  coordinate  displacements, 
as  defined  by  Equations  (4.4)  and  (5.  12),  let 


(“al} 

I  r“a2>  [ 

.  C“a3>„ 


(5.24) 


Consistent  with  Equation  (5.  10),  let  the  consolidated  vector  of  system 
coordinate  displacements  be 


<“0  ■ 


(“a.) 

{‘a2> 

(aa3> 


< 


(5.25) 


Corresponding  to  Equations  (5.15),  (5.  24)  and  (5.25),  the  transformation 
matrix  ItJ  is  defined  by 


where 


{»b)=  [Tb]K> 


(5.26) 


(5.27) 


From  Equations  (5.  10)  and  (5.25),  it  is  noted  that  the  coordinate  identifica¬ 
tion  vector  for  fluid  component  b  is  obtained  as  a  combination  of  the  shell 
coordinate  identification  vectors  defined  in  Equation  (5.  11),  for 

a  =  al,  a2  and  a3,  as  follows: 


<Cal) 


<Cb)  *  <  (Ca2>  } 


(5.28) 


5.  9 


Table  1.  Shifted  Chebyshev  Polynomial  Coefficients  for  Use  i 
Constructing  Matrices  and  jj3j 
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6.  LAUNCH  VEHICLE  STIFFNESS  AND 
MASS  MATRIX  SYNTHESIS 

In  order  to  construct  the  total  vehicle  stiffness  and  mass  matrices 
referenced  to  a  common  coordinate  system,  the  individual  component 
stiffness  and  mass  matrices  are  expressed  in  terms  of  the  system  coor¬ 
dinates  developed  in  Section  5,2.  For  each  shell  and  fluid  component  e, 

4 

the  matrices,  transformed  to  system  coordinates,  become 

[Ke]  -  [Te]Tfc]  [Te] 

and  (6.1) 

M =  NT  H  [Te] 

where  T  is  the  transformation  matrix  defined  in  Equation  (5.  13)  for 

L  eJ  r  -i 

shell  components  and  in  Equation  (5.  26)  for  fluid  components.  K  and 

r—  i  l  e  j 

Ml  are  the  stiffness  and  mass  matrices  related  to  the  local  coordinate 
system.  The  stiffness  and  mass  matrices  for  the  spring  mass  components 
do  not  undergo  the  transformation  Equation  (6.  1),  since  they  already 
exist  as  input  data  in  terms  of  the  system  coordinate  displacements. 

The  total  system  stiffness  and  mass  matrices  are  synthesized  by 
expanding  each  of  the  component  matrices  Equation  (6.  1)  into  an  enlarged 
matrix  which  is  of  the  same  order  as  the  total  system  matrix  and  which 
is  related  to  the  set  of  system  coordinates  for  the  total  vehicle,  as  shown 
in  Figure  11.  This  is  accomplished  by  superimposing  each  of  the  com¬ 
ponent  matrices  after  an  additional  transformation,  as  described  below. 

As  indicated  in  Figures  8,  9,  10  and  1 1  and  discussed  in  Section  5.  2, 
each  system  displacement  and  rotation  coordinate  is  defined  by  an  identi¬ 
fication  number  C(C  =  1,2,...,  N^)  arbitrarily  specified  to  suit  the 
convenience  of  the  analyst.  It  is  desirable,  however,  to  reserve  the 

identification  numbers  (N  ),  (N  -  1),  .  .  .  ,  (N  -  N  +1)  for  the  N 

c  c  c  o  o 

system  coordinates  to  be  used  for  rigid  supports  to  preserve  a  one-to-one 
correspondence  in  the  digital  program  output  between  the  identification 
number  s  and  the  coordinate  row-column  location  in  the  total  system  matrices. 


The  arrangement  of  the  coordinates  in  the  total  system  displacement 
vector  {«}  of  order  is  made  to  coincide  with  the  coordinate  identi¬ 
fication  numbers  C.  Thus,  the  transformation  matrix  M  relating 
the  total  system  displacements  to  the  component  e  system  displace¬ 
ments  is  defined  as 


(“e)  =  [AJ{“}  '  (6'21 

The  elements  of  are  determined  to  be  zero  or  unity  from 


the  relations 


6  =1  when  s  =  C 

r  s  r 


=  0  when  s  ^  C 


(6.3) 


where  C  is  the  r-th  element  of  the  coordinate  identification  vector  f  C  > 

r  r  i  V  e/ 

for  component  e  [see  Equations  (5.11)  and  5.  28)  « 

The  total  stiffness  and  mass  matrices  are  obtained  using  transfor¬ 
mation  matrices  M  in  the  following  summation  equations 


C-Wv  -i  H’  MM  *.£  WT HW 


”C-”o'  WT  WM  *j'WTWW 


M 

+  l. 


Swtww 


(6.5) 


in  which  N^,  Ng,  and  N^.  are  the  total  number  of  system  coordinates, 
shell  components,  fluid  components  and  spring -mass  components,  re¬ 
spectively.  The  matrices  JkJ  ,  and  are,  respectively,  the 

stiffness  and  mass  matrices  for  the  shell  components  and  the  mass  matrix 
for  the  fluid  components,  as  defined  by  Equation  (6.  1). 


The  superposition  technique  assures  displacement  compatibility 
and  force  equilibrium  at  the  joints  between  components.  Displacement 
boundary  conditions  are  imposed  on  the  total  stiffness  and  mass  matrix 
by  removing  appropriate  rows  and  columns  of  coefficients  corresponding 
to  points  on  the  vehicle  and  its  support  which  are  rigidly  restrained  from 
motion.  Due  to  storage  limitations  in  the  digital  program,  the  order  of 
the  resulting  total  stiffness  and  mass  matrix  must  not  exceed  80. 


Additional  limitations  are  placed  on  the  size  of  the  component 
matrices  utilized  in  Equations  (6.4)  and  (6.5).  The  shell  component 
stiffness  and  mass  matrices,  and  must  have  an  order  no 

larger  than  22.  For  the  fluid  mass  matrix,  J  ,  the  order  must  not 
exceed  35,  and  for  the  spring -mass  component  stiffness  and  mass 
matrices,  JkJ  and  |\tcJ  ,  the  order  must  not  exceed  10. 


Figure  11.  Vehicle  System  Coordinates 
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7.  DYNAMIC  RESPONSE  EQUATIONS 


7.  1  NATURAL  FREQUENCY  EQUATIONS 

The  total  stiffness  and  mass  matrices  which  are  derived  in 
Section  6  are  used  for  computation  of  the  natural  frequencies  and  mode 
shapes  from  the  eigenvalue  equation 


[k]  (a)  -  p2  [m]  {a}  =  0 


(7.  1) 


in  which  p  is  the  circular  frequency  of  the  launch  vehicle  and  {*)  is  the 
modal  vector  whose  elements  are  the  longitudinal,  radial  and  rotational 
system  coordinate  displacements  defined  in  Section  5  and  illustration  in 
Figure  11.  This  equation  is  solved  to  obtain  the  natural  frequencies  p^ 
and  the  mode  shapes  for  all  modes,  t,  which  are  arranged  in  a  square 
modal  matrix  of  order  N^  -  Nq.  Each  column,  {“*)•  of  ,  is  the 

mode  t  displacement  vector  with  system  coordinate  elements,  whereas 
each  row,  |jagJ  ,  °f  [aj  ,  is  the  system  coordinate  s  displacement  vector 
with  natural  mode  elements. 


7.2  STEADY -STATE  RESPONSE  EQUATIONS 

The  steady-state  response  due  to  simple  harmonic  loads  of  fre¬ 
quency  co  is  determined  using  a  standard  modal  technique.  The  elements 
of  the  load  vector  (p)  of  order  (Nc  -  Nq)  represent  axisymmetric  forces 
(longitudinal  and  radial)  or  moments,  depending  on  whether  the  associated 
coordinate  is  a  displacement  or  a  rotation.  The  displacement  response 
{r}  at  coordinates  s(s  =  1,  2,  ....  (Nc  -  Nq)^  on  the  launch  vehicle  is 
expressed  as  the  linear  superposition  of  the  individual  modal  responses 
based  on  an  assumed  modal  damping  factor  which  is  the  ratio  of  the 
actual  damping  to  the  critical  damping  for  each  mode  and  has  the  form^ 


In  these  equations,  6  =  tt  when  p^_  =  0,  and  6^  =  0  when  0,  to  =  0. 

(r)  and  {b)  represent,  respectively,  the  vectors  of  steady-state  dis¬ 
placement  amplitude  and  the  phase  angle  by  which  the  forcing  function 
leads  (+)  or  lags  (-)  the  response.  The  velocity  and  acceleration 

<*>  responses  are  obtained  from  the  relations 


co  <  R  sin 


(wt '  (r +  6))} 


(r)  =  -w2  (r  sin  (wt  -  6)} 


(7.5) 


The  internal  forces  (or  moments)  <Sa>  acting  at  each  point  along 
the  vehicle  on  each  shell  component  a  are  obtained  from  the  equation 


(s  )  =  (s  sin  (wt  -  6)} 

V  a/  ((U+V)  xl  ) 


(7.6) 
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where 


02)((U+V>xl)  =  ([Ka]  H  (R  sin  ?>)2  +  (Kj  [aJ  (R  cos  5})' 


[(U+  V)  x  l) 


=  /  tan 


-ifjgal  M(R  sin6V 

usJ  [Aa](fcos^) 
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and  ^5}  represent,  respectively,  the  amplitude  and  phase  angle  of 
the  internal  forces.  |_^SJ  row  s  (s  =  2,  •  *  “  ,  (U  +  V)a)  of  the  shell 

element  stiffness  matrix  Kl-  00  and  CO  are  obtained  from 
Equation  (7 . 3). 


8.  DIGITAL  PROGRAM  ARRANGEMENT 


The  analytic  model  developed  in  the  previous  sections  is  used  as 
the  basis  for  a  digital  computer  program  to  determine  the  vibration 
characteristics  and  steady-state  response  of  a  launch  vehicle  subjected 
to  axisymmetric  sinusoidal  loads.  The  program  is  written  in  Fortran  IV 
language  for  use  on  an  IBM  7094  computer  having  32K  magnetic  core 
storage  locations.  The  functional  operations  and  the  overall  program 
arrangement  are  illustrated  in  Figure  12.  Each  of  the  operations  en¬ 
closed  by  a  block  represents  an  independent  link  in  the  computer  program. 
This  feature  facilitates  the  task  of  making  future  modifications  or  expan¬ 
sions  of  the  program. 

The  input  data  necessary  to  inititate  the  program  sequence  is  dis¬ 
cussed  in  detail  in  Section  9.  After  the  data  has  been  processed,  the 
program  proceeds  to  develop  the  stiffness  and  mass  matrices  for  each 
of  the  shell  components.  The  functional  operations  required  to  perform 
these  computations  and  the  form  of  the  equations  used  in  the  digital  pro¬ 
gram  are  provided  in  Appendix  D.  Consistent  with  the  matrix  formulation 

of  the  basic  equations,  the  numerical  integration  is  performed  using  a 
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sixteen  point  Gaussian  weighted  matrix  integration  scheme. 

In  a  similar  fashion,  the  fluid  mass  matrix  for  each  fluid  compo¬ 
nent  is  constructed.  The  functional  flow  diagram  and  the  basic  equations 
required  to  describe  the  three  fluid  tank  configurations  considered  in  the 
present  analytical  model  are  presented  in  detail  in  Appendix  E.  Unlike 
the  formulation  for  the  shell  stiffness  and  mass  matrices,  the  equations 
for  the  fluid  mass  matrix  involve  a  double  integration.  For  this  compu¬ 
tation,  a  double  Lagrangian  weighted  matrix  integration  scheme  was 
found  most  suitable.  This  technique  employs  two  11 -point  Lagrangian 
weighting  matrices  in  sequence  to  provide  a  22  point  approximation. 

The  component  matrix  construction  is  concluded  with  the  setting  up 
of  the  stiffness  and  mass  matrices,  provided  as  input  data,  of  the  spring - 
mass  components. 


The  shell,  fluid  and  spring-mass  component  stiffness  and  mass 
matrices  are  then  synthesized  into  a  total  vehicle  system  stiffness  and 
mass  matrix,  according  to  the  steps  presented  in  Section  6.  Utilizing 
these  matrices,  the  natural  frequency  equation,  Equation  (7.  1),  is  for¬ 
mulated  and  subsequently  solved  using  a  standard  digital  eigenvalue 
routine  which  solves  matrices  up  to  order  80, 

The  program  user  now  has  the  option  of  1)  continue  the  analysis 
and  move  directly  to  the  computation  of  the  steady-state  response,  or 
2)  to  temporarily  stop  the  solution  after  the  free  vibrations  stage  for 
the  purpose  of  examining  the  output  before  proceeding  with  the  compu¬ 
tation  of  the  steady-state  response.  This  option  enables  the  user  to 
examine  the  results  of  the  computation  for  the  natural  frequencies  and 
mode  shapes  before  determining  the  input  for  the  modal  damping  and 
for  the  frequencies  at  which  the  launch  vehicle  is  forced. 

Checks  on  the  accuracy  and  consistency  of  the  various  computations 
are  performed  throughout  the  program.  The  shell  component  stiffness 
matrices  are  subjected  to  longitudinal  rigid  body  displacements  to  estab¬ 
lish  automatically  that  equilibrium  is  satisfied.  For  each  shell  compo¬ 
nent,  the  "equilibrium  check"  appears  as  output  in  the  following  form 


Equilibrium 

Check 


/L  Longitudinal  Rigid  Body  Displacement  Forces 

U 

£ 

i=l 


k.. 

li 


(8.  1) 


where  k^  are  the  diagonal  components  of  the  shell  component  stiffness 
matrix  |kJ  associated  with  the  system  coordinates  and  U  is  defined 
in  Section  5.2  as  the  total  number  of  longitudinal  system  coordinates 
associated  with  component  a.  If  equilibrium  is  satisfied,  Equation  (8.  1) 
will  be  equal  to  zero. 

In  a  similar  manner,  the  mass  matrices  for  the  shell  and  fluid 
components  are  subjected  to  unit  accelerations  to  verify  that  the  rigid 
body  mass  characteristics  have  been  accurately  represented.  The 
"mass  check"  also  appears  as  output  and  has  the  form 

M. 

input 
M 


Mass  Check  = 


computed 


(8.2) 


where  M^nput  is  th e  total  mass  of  the  element,  provided  as  input,  and 
^"computed  body  mass  computed  in  the  program.  The  ele  - 

ments  of  the  mass  matrix  [mJ  are  corrected  by  the  "mass  check" 
factor  to  provide  the  correct  total  mass  representation  for  the  structural 
component  e,  since  ^computecp  general,  will  not  agree  with  Mqnput 
due  to  the  physical  nature  of  the  problem. 

An  additional  check  on  the  accuracy  of  the  program  is  provided  in 
the  form  of  a  FtJ  "inverse  check."  In  developing  the  transforma¬ 

tion,  Equation  (5.  15),  it  is  necessary  to  take  the  inverse  of  £*T J  An 
indication  of  the  conditioning  of  this  matrix  is  furnished  by  the  product 


Inverse  Check  = 


X]  [’ 


“  -1 
a 


(8.3) 


which  appears  as  output.  The  deviation  of  this  product  from  a  unit  diago¬ 
nal  matrix  provides  an  estimate  of  the  accuracy  of  the  computation  for 
£tJ  .  The  incorporation  of  the  three  checks  discussed  above  is  a  useful 
aid  in  assuring  the  reliability  of  the  digital  solution. 

An  additional  feature  which  is  incorporated  in  the  program  is  the 
capability  to  stack  cases,  i.  e.,  to  solve  numerous  related  problems  in 
sequence  by  simply  stacking  the  data  input  for  each  of  the  cases.  The 
stacking  capability  can  be  utilized  for 

1)  the  solution  of  the  complete  problem,  which  may  include 
the  steady-state  response  or 

2)  the  solution  of  the  steady-state  response  using  natural 
mode  data  previously  stored  on  tape. 

This  feature  provides  for  an  efficient  use  of  machine  time  by  eliminating 
the  need  for  reloading  the  program  deck  for  successive  cases. 


CONSTRUCT: 

TOTAL 

[K]  AND  [m] 
MATRICES  FROM 
COMPONENT  MATRICES 


COMPUTE: 
FREQUENCIES, 
MODE  SHAPES, 
VELOCITIES  AND 
ACCELERATIONS 


CONSTRUCT: 


N 

FOR  N 


AND 


M 


M 


spring- 


mass  COMPONENTS 


COMPUTE: 
STEADY-STATE 
RESPONSE  DISPLACEMENTS, 
VELOCITIES,  ACCELERA¬ 
TIONS,  AND  COMPONENT 
INTERNAL  FORCES 


Figure  12.  Digital  Program  Arrangement 


9.  DATA  SETUP 


9 . 1  VEHICLE  SUBDIVISION 

In  order  to  prepare  the  input  data,  it  is  first  necessary  to  subdivide 
the  launch  vehicle  into  a  consistent  set  of  axisymmetric  shell  components 
a,  fluid  components  b,  and  spring- mass  components  c.  Size  limitations 
of  the  program  require  that 

1)  The  total  number  of  shell  components  shall  not  exceed 
forty  (40), 

2)  The  total  number  of  fluid  components  shall  not  exceed 
six  ( 6) , 

3)  The  total  number  of  spring- mass  components  shall  not 
exceed  thirty  (30),  and 

4)  The  order  of  the  spring-mass  stiffness  and  mass  matrices 
must  not  exceed  ten  (10). 

Subsequent  to  the  vehicle  subdivision,  the  location  of  the  system 
coordinate  displacements  must  be  determined.  The  internal  operations 
of  the  digital  program  require  that 

1)  Longitudinal,  radial  and  rotational  coordinate  displacements 
must  be  placed  at  each  junction  of  two  or  more  shell  elements 

2)  Longitudinal  and  radial  coordinate  displacements  are  placed 
at  each  point  lying  on  the  vehicle  longitudinal  axes. 

These  locations  are  defined  as  the  terminal  points  of  each  shell 
element.  Additional  longitudinal  and  radial  coordinate  displacements  may 
be  placed,  as  desired,  at  intermediate  points  uniformly  spaced  between 
the  terminal  points.  Longitudinal,  radial  and  rotational  coordinates  may 
also  be  used  to  describe  the  motions  of  the  spring-mass  elements. 

The  coordinate  displacements  are  then  identified  by  a  consecutively 
numbered  sequence.  The  arrangement  of  these  numbers  on  the  vehicle  is 
arbitrary  and  is  left  up  to  the  discretion  of  the  analyst.  However,  the 
radial  coordinate  points  lying  on  the  vehicle  axis  are  not  to  be  identified 
with  a  number.  This  is  necessary  because  the  program  must  assume  these 
displacements  equal  to  zero,  as  is  required  by  the  axisymmetry  of  the  ve¬ 
hicle,  and  thus  they  do  not  contribute  to  the  degrees -of-freedom  of  the 
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system.  These  are  defined  as  "unnumbered”  coordinates  and  are  utilized 
only  in  the  preparation  of  the  data  input  sheets  for  determining  the  value 
of  U  and  V  (see  Figure  11). 

In  general,  for  a  supported  structure,  some  of  the  "numbered"  co¬ 
ordinates  may  be  assumed  fixed  or  restrained  from  motion.  The  total 
number  of  "numbered"  coordinates  which  are  not  fixed  must  not  exceed 
eighty  (80).  In  addition,  the  numbering  sequence  should  be  arranged  so 
that  the  fixed  coordinates  are  numbered  last  to  insure  consistency  in  the 
program  output  identification  with  the  coordinate  identification  numbers. 

The  total  number  of  coordinates  for  a  shell  component  may  not  exceed  22, 
of  which  in  general  the  number  of  longitudinal  displacement  coordinates  may 
not  exceed  11,  and  the  number  of  lateral  displacement  coordinates  may  not 
exceed  10  (see  Section  5.3).  For  the  case  of  three  shell  components  sur¬ 
rounding  a  fluid  component,  the  total  number  of  coordinates  for  the  com¬ 
posite  tank  structure  may  not  exceed  29. 

9.2  INPUT  DESCRIPTION 


1)  Heading 
HHEAD 


System  Input  Data 


is  one  line  of  BCD  characters  which  will 
be  printed  as  the  title  of  the  printed  out¬ 
put.  The  number  of  BCD  words  (6  char¬ 
acters  per  word)  must  not  exceed  11. 


2)  Input  Parameters 

N~  is  the  total  number  of  system  coordinates 

which  include  the  fixed  coordinates  N^  . 
N_  -  N0  must  not  exceed  80. 

\-J 

Nq  is  the  total  number  of  shell  components  . 

Nc  cannot  be  zero  and  must  not  exceed  40. 
o 

N  is  the  total  number  of  fluid  components. 

F 

N^  must  not  exceed  6. 

N^.  is  the  total  number  of  spring-mass  com¬ 

ponents.  N^.  must  not  exceed  30. 
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Nq  is  the  total  number  of  fixed  coordinates. 

Nc  -  N q  must  not  exceed  80. 

3)  Applied  Loads  and  Forcing  Frequencies  (Section  7.2) 

Nt  is  the  total  number  of  discrete  applied 

L 

loads.  must  not  exceed  80. 

C.  is  the  applied  load  coordinates  C  ,  C  , 

X  L  Ci 

Cy  •  •  •  ,  ClSlL- 

P.  is  the  discrete  applied  loads  P  ,  P  ,  P  , 

l  1  u  j 

•••>  PNL- 

N^.  is  the  number  of  sets  of  forcing  function 

frequencies  . 

L,  Af.,  rm  f  is  the  frequency  of  the  forcing  function 

in  cycles  per  second  (to  =  Zirf ) .  This  program 

will  compute  the  steady- state  response  for  the 

frequencies  f . ,  f.  +  Af . ,  f.  +  2  Af . ,  *  *  *  , 
1111  i 

f.  +(m.  -  1)  Af.,  is  1,  2,  •••,  Nw.  ' 


4)  Modal  Damping  Factor 


is  the  number  of  input  r\ .  Program  will 
generate  a  complete  table  of  r|  by  setting 

T1(Net  +  1)’  r|(NET+2)’  ‘  ‘  ’  T,(Nc-N0)  equal 

^1  NET 

is  the  ratio  of  the  assumed  damping  to  the 
critical  damping  in  mode  k 


^2*  ^3'  ’ 


5)  Ratio  of  Acceleration 


is  the  ratio  of  the  vehicle  acceleration  to 
the  acceleration  of  gravity. 
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6)  Steady-State  Response  Option 

S  is  a  fixed  point  word  which  controls  the 

option  of  computing  steady-state  response. 

S  0  indicates  that  the  computation  of  the  steady- 

state  response  is  not  included.  The  nec¬ 
essary  data  for  the  steady- state  computa¬ 
tion  is  saved  on  Tape  1  and  2. 

S  =  0  indicates  that  the  computation  of  the  steady- 

state  response  is  included. 

S  0  is  the  option  to  compute  the  steady- state 

response  only.  The  necessary  data  should 
be  available  on  Tape  1  and  2.  In  addition 
the  following  System  Input  Data  must  be 
provided: 

S  =  1,  Heading  [item  1)],  Input  Param¬ 
eters  [Item  2)],  and  Modal  Damping 
Factors  [item  4)]. 

S  =  2,  Heading  [Item  1)],  Input  Param¬ 
eters  [Item  2)],  Applied  Loads  and 
Forcing  Frequencies  [Item  3)], 

Modal  Damping  Factors  [item  4)], 
and  the  option  word  opt^. 


7)  Print  Options 


opt^  is  an  option  word  which  controls  the  output 

of  stiffness  matrix  and  mass  matrix  of  the 
shell  and  the  fluid  components. 

opt^  =  1,  print  the  component  matrices 

opt  =  0,  suppress  the  printing  of  com¬ 
ponent  matrices. 

opt^  is  an  option  word  which  controls  the 

printing  of  total  stiffness  matrix  and  total 
mass  matrix. 

opt^  =  1,  print  the  total  stiffness  and  mass 
matrices. 
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opt^  =  0,  suppress  the  printing  of  total 
stiffness  and  mass  matrices. 

is  an  option  word  which  sets  the  rigid 
body  frequency  to  zero  for  computing  the 
response . 

opt^  =  1,  set  the  first  frequency  to  zero. 

opt^  =  0,  do  not  set  the  first  frequency  to 
zero. 

is  an  option  word  which  controls  the  com¬ 
putation  and  printing  of  forces  for  the 
steady-state  response. 

opt4  =  1,  compute  and  print  the  forces. 

opt^  =  0,  do  not  compute  and  print  the 
forces. 

is  the  number  of  frequencies,  mode  shapes, 
velocities  and  accelerations  that  will  be 
printed  as  the  final  output. 


8)  Polynomial  Matrices  (Section  5 . 3) 


is  the  total  number  of  polynomial  matrices, 
is  the  number  of  rows  of  polynomial  matrix 
M  Ukxl  1  •  <0k  =  "> 

is  the  number  of  rows  of  polynomial  matrix 

M  v.xir  = 


is  U^_  x  11  polynomial  matrix. 


is  x  11  polynomial  matrix, 
k  =  1,  2,  3,  •  • 


The  input  sequence  of  the  polynomial  ma¬ 
trices  establishes  the  identification  number 
k  which  is  referred  by  the  shell  components. 
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The  subscript  k  is  used  as  the  polynomial 
matrix  identification  number  by  the  shell 
components . 


Shell  Component  Input  Data 


1)  I.D.  Number 


2)  Coordinates 
U,  V 
U,  V 


is  the  identification  number  for  shell  com¬ 
ponent  a  where  0  a  Ng 

4*  a  indicates  a  conical  shell  component 
-  a  indicates  an  ellispoidal  shell  component 


are  the  total  number  of  system  coordinates, 

are  the  total  number  of  local  coordinates 
U,  V  must  not  exceed  1 1 . 


3)  Coordinate  I.D.  Vector  (Figures  8,  9,  and  10) 

(ID).  is  the  identification  vector  which  is  used 

to  position  the  elements  for  building  total 
stiffness  and  mass  matrices.  The  length 
of  the  vector  must  be  equal  to  U  +  V  and 
the  number  must  not  be  greater  than  N^,. 

4)  Polynomial  Matrix  Identification  Number 

k  is  the  polynomial  matrix  identification 

number  which  refers  to  polynomial  matri¬ 


ces 


|aJ  and  £bJ  k  in  the  system  input  data. 


5)  Shell  Geometric  Data  (Figures  4  and  5) 

c|>o  is  the  meridional  angle  for  conical  shell 

and  is  the  edge  meridional  angle  for  ellip¬ 
soidal  shell.  cj>  is  input  in  degrees. 

L  is  the  height  of  conical  shell 

+  L  indicates  converging  upward 
-  D  indicates  converging  downward 
L  =  0  for  ellipsoidal  shell  input 
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is  the  lower  radius  of  conical  shell 
=  0  for  ellipsoidal  shell  input 

b  is  the  height  of  ellipsoidal  shell 

+  b  indicates  convex  upward 
-  b  indicates  convex  downward 
b  =  0  for  conical  shell  input 

a  is  the  radius  of  the  base  of  ellipsoidal 

shell,  a  =  0  for  conical  shell  input 

6)  Orthotropic  Shell  Constants  and  Thickness  (Equation  3.5) 


<CU>p 

<C12>p 

<C22»p 

|C33>p 

<C34>p 

<C44>p 


are  orthotropic  shell  constants  at  two 
points  £  =  0,  1  which  are  represented  by 
p  =  1,  2,  respectively. 

are  orthotropic  shell  constants  at  four 
points  £  =  0,  1/3,  2/3,  1  which  are 
represented  by  p  =  1,  2,  3,  4,  respectively. 

are  shell  thickness  at  two  points  £  =  0,  1 
which  are  represented  by  p  =  1,  2, 
respectively. 


7)  Mass  Density  and  Total  Mass 


is  the  mass  density  of  the  shell  component. 


is  the  total  mass  of  the  shell  component. 

M  ^  0,  the  ratio  of  the  total  mass  M  to 
a  __  a 

the  computed  mass  M  will  be  used  as  the 

a 

scaling  factor  for  the  mass  matrix.  When 
opt^  =  1,  the  scaling  factor  will  be  printed 
as  the  mass  check  of  the  mass  matrix. 

M^  =  0,  no  scaling  factor  will  be  used  for 
the  mass  matrix.  When  opt^  =  1,  the 
computed  mass  will  be  printed. 
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8)  Initial  Stress  Data  (Figures  (B.  1),  (B.  2)  and  )B.  3)) 


H. 

x 


w. 

i 


W 


is  the  depth  of  interior  fluid. 

is  the  weight  density  of  interior  fluid. 

is  the  uniform  interior  pressure. 

is  the  depth  of  exterior  fluid. 

is  the  weight  density  of  exterior  fluid. 

is  the  uniform  exterior  pressure 

is  the  reactive  force  at  upper  edge  of 
conical  shell. 

+  W  produces  tensile  stresses. 

-  W  produces  compressive  stresses. 
W  =  0  for  ellipsoidal  shells. 


Fluid  Component  Input  Data 

1)  I.D.  Number 

b  is  the  identification  number  for  fluid 

component  b  where  0  b  Zf  N  . 

-T 

2)  Associated  Shell  Components  (Figure  7) 

a  ,  a  ,  and  are  the  identification  numbers  of  the 

1  £ 

a^  associated  shell  components 


3)  Fluid  Data  (Figure  C.  1) 

H  is  the  depth  of  fluid  component. 

y  is  the  mass  density  of  fluid  component. 

M  is  the  mass  of  fluid  component.  M  =£  0, 

the  ratio  of  the  total  mass  M  to  the  com¬ 
puted  mass  M  will  be  used  as  the  scaling 
factor  for  the  mass  matrix.  When  opt ^  =  1 , 
the  scaling  factor  will  be  printed  as  the 
mass  check  of  the  mass  matrix. 
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M  =  0,  no  scaling  factor  will  be  used  for 
the  mass  matrix.  When  opt^  =  1,  the 
computed  mass  will  be  printed. 


Spring-Mass  Component  Input  Data 

1)  LD.  Number 

c  is  the  identification  number  of  spring- 

mass  component  c  where  0  -*=  c  ^  N^.. 


2)  Stiffness  and  Mass  Matrices 


Me 

Me 


is  the  order  of  the  spring- mass  component 
n  must  not  exceed  10. 

is  n  x  n  stiffness  matrix  of  spring- mass 
component. 

is  n  x  n  mass  matrix  of  spring -mass 
component . 


3)  Coordinate  I.D.  Vector 


IDC 

i 


is  the  identification  vector  which  is  used 
to  position  the  elements  for  building  the 
total  stiffness  and  mass  matrices.  The 
length  of  the  vector  must  be  equal  to  n 
and  the  number  must  not  be  greater  than 


N 


C' 


9.3  SAMPLE  INPUT  DATA  SHEETS 


Sample  input  data  sheets  are  included  to  illustrate  the  input  format. 
The  data  must  be  arranged  in  the  order  shown,  that  is, 

1)  System  input  data 

2)  Shell  component  input  data 

3)  Fluid  component  input  data 

4)  Spring -mass  component  input  data 
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THE  SYSTEM  INPUT  DATA 


SPACE  TECHNOLOGY  LABORATORIES,  INC. 

OATE _  COMPUTATION  AND  DATA  REDUCTION  CENTER  PAOE _ OF _ 


9.10 


9.12 


THE  SHELL  COMPONENT  INPUT  DATA 
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DATE 


SPACE  TECHNOLOGY  LABORATORIES,  INC. 
COMPUTATION  AND  DATA  REDUCTION  CENTER 


PAGE 


THE  FLUID  COMPONENT  INPUT  DATA 


9-15 


THE  SPRING-MASS  COMPONENT  INPUT  DATA 
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10.  NOTATION 


Row  matrix  of  order  j 

Column  matrix  of  order  i 

Rectangular  matrix  with  i  rows  and  j 
columns 


a,  b,  c  Identification  number  for  the  shell,  fluid  and 
spring-mass  components  respectively 

a,  b  Semimajor  and  semiminor  axes,  respectively, 
of  an  ellipsoidal  bulkhead 

a,  T>  Edge  radius  and  height,  respectively,  of  an 
ellipsoidal  bulkhead 


al,  aZ,  a3 

M(Uxll)  =  phcn] 
[B](Vxll)  [bi  n] 
(Ce> 

C 

r 

Cll«>-  C12®-  C22<» 

c33(4),  c34(4),  c44(|) 

ds 


Identification  numbers  for  the  shell  components 
which  enclose  a  fluid  component 

Polynomial  matrix  associated  with  u^d) 

Polynomial  matrix  associated  with  v^  (|) 

Coordinate  identification  vector  for 
component  e 

r*k  component  of 

Orthotropic  stress- strain  coefficients 

Orthotropic  moment- curvature  coefficients 
Differential  meridional  distance  along  shell 

Derivative  with  respect  to  meridional 
distance  s 


Dj,  D 


Z 


Constants  used  to  determine  the  rotation 


vector 


[D,l  ,  [DJI  Matrices  used  in  the  definition  of  the 

^  J(Zxll)  -'(Zxll)  rotation  vector  (p\ 
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e  General  identification  number  for  vehicle  com¬ 
ponents  which  may  stand  for  a  or  b 

g  Ratio  of  vehicle  acceleration  to  acceleration 
of  gravity 


H  Total  fluid  level  measured  positive  upward 
from  the  base  of  a2 


Hf,  H^,  Fluid  levels  associated  with  tank  shell 

sections  al,  a2,  a3 


[K](N  xN  ) 
'  C  c' 

J  ((U+V)x(U+V)) 
j  ((U+Y)x(U+  V)) 


Total  launch  vehicle  stiffness  matrix 


Stiffness  matrix  for  shell  component  a 
associated  with  the  system  coordinates 


Stiffness  matrix  for  shell  component  a 
associated  with  the  local  coordinates 


Stiffness  matrix  of  spring- mass  component  c 
associated  with  the  system  coordinates 


Row  s  (s  =  1,  2,  •  •  •  (U  +  V)J  of  [kj 


K,  Shell  meridional  curvature 

0 


Ke 

L 


Shell  hoop  curvature 

Length  of  conical  shell 

Generalized  mass  for  mode  t 

Total  mass  matrix  for  the  launch  vehicle 

Mass  matrix  for  shell  component  a 
associated  with  the  system  coordinates 

Mass  matrix  for  shell  component  a 
associated  with  the  local  coordinates 

Mass  matrix  for  fluid  component  b 
associated  with  the  system  coordinates 
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py 

N 

NC 


N 


O 


Ns 


p 

Pt 

(p>  ((Nc-N0)Kl) 


Mass  matrix  for  fluid  component  b 
associated  with  the  local  coordinates 

Mass  matrix  for  spring-mass  component  c 
associated  with  the  system  coordinates 

Total  number  of  system  coordinates  used  in 
the  vehicle  model 

Total  number  of  fluid  components  used  in 
the  vehicle  model 

Total  number  of  spring-mass  components 
used  in  the  vehicle  model 

Total  number  of  fixed  coordinates 

Total  number  of  shell  components  used  in 
the  vehicle  model 

Initial  meridional  stress 

Circular  frequency  of  the  vehicle 

Circular  frequency  of  the  vehicle  for  mode  t 

Applied  load  vector 

Generalized  force  acting  on  mode  t 

Radial  distance  from  the  vehicle  longitudinal 
axis  to  each  point  on  the  shell 

Radial  distance  from  the  vehicle  longitudinal 
axis  to  each  point  in  the  fluid 

Meridional  radius  of  curvature  of  the  shell 

Hoop  radius  of  curvature  of  the  shell 

Displacement  response  vector 

Displacement  response  amplitude  vector 
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Vector  of  internal  forces  or  moments  acting 
at  each  point  of  shell  component  a 


[T a)  ((U+V)x(U+V)j 

£  *3  (WxW) 
(( U+  V)x  l)  ’  ({ U+  V)x  l) 

^(S).  v(e) 

("(x)}(Wxl)={v(x)} 

and 

{v(x,  r)}(Wxl)  ={vm(x,  r)} 

U,  V 

U,V 

M(Uxll)’  M(Vxll) 


Time 


Transformation  matrix  which  relates  local 
to  system  coordinates  in  shell  component  a 

Transformation  matrix  which  relates  local 
to  system  coordinates  in  fluid  component  b 

Longitudinal  and  radial  displacement  vectors, 
respectively,  for  shell  components 

Generalized  longitudinal  and  radial  displace¬ 
ments,  respectively,  for  shell  components 

Longitudinal  and  radial  fluid  displacements, 
respectively 

Total  number  of  longitudinal  and  radial 
system  coordinates,  respectively,  associated 
with  each  shell  component 

Total  number  of  longitudinal  and  radial  local 
coordinates,  respectively,  associated  with 
each  shell  component 

Constant  matrices  used  in  the  definition  of 
the  system  coordinates 


w.,  w  Fluid  weight  densities  interior  and  exterior, 

respectively,  to  each  shell  component 


w.,  w  Effective  fluid  weight  densities  equal  to, 

respectively,  gw^  and  gw^ 

W  The  sum  of  U  +  V  for  the  three  shells 
surrounding  a  fluid  component 
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x  Longitudinal  axis  of  the  launch  vehicle 


z Structural  impedance  for  mode  t 


Modal  vector  whose  components  are  the 
longitudinal,  radial  and  rotational  system 
coordinate  displacements  for  shell  com¬ 
ponent  a 


^  a^((U+V)xl) 


(Uxl) 


ih)  - 

v  J  (Vxl) 


Consolidated  vector  of  local  coordinates 


Generalized  coordinates  in  the  longitudinal 
and  radial  directions,  respectively 


h 


(Nc-N0)k(nc-no) 


Modal  matrix  with  each  column  representing 
a  vehicle  mode  shape 


Each  row  of  £aj  which  is  a  system  coordinate 
"s"  displacement  vector  with  natural  mode 
elements 


(NC“N0)xl 


Each  column  of  [a]  which  is  a  mode  Mtn 
displacement  vector  with  system  coordinate 
elements 
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a* 


Mass  densities  for  shell  and  fluid  components 
respectively 


Vector  of  steady- state  displacement  phase 
angles 


((U+V)xl) 


Vector  of  phase  angles  for  the  internal 
forces  of  each  shell  component 


M  ((U+V)xNc)  "  [5rs) 


Transformation  matrix  relating  total  system 
displacements  to  component  3  system 
di  splac  ement  s 


€ 


€e 


Meridional  and  hoop  strains,  respectively, 
for  the  shell  components 


'H 


t 


Ratio  of  actual  damping  to  the  critical  damping 
for  each  mode  t 


Nondimensional  variable  describing  location 
on  each  shell  component 
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p  Meridional  rotation  of  the  shell  components 

<p  Meridional  angle 

co  Frequency  of  the  forcing  function  in  radians 
per  second 
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APPENDIX  A 

SHELL  GEOMETRY  AND  ENERGY  EXPRESSIONS 


The  additional  potential  energy  of  a  shell  of  revolution  due  to  axisym- 
metric  deformations  is  given  in  the  form 


V 


=  1  f^(: 


N  €  „  +  N„cq  +  M  K  +  M.KL  +  N 


0  0 


He'e 


0  0 


0  0 


T  °  ^ 

V p ) 


ds 


(A.  1) 


in  which  the  last  term  represents  the  work  done  by  the  initial  meridional 

stress,  N  °.  ^  The  initial  hoop  stress  does  not  make  a  similar  contribution 

0 

to  the  potential  energy  since  there  is  zero  rotation  in  the  hoop  direction. 

g 

In  the  notation  of  Flugge,  the  strains  €q),  curvatures  (K^,  Kg)and  the 

meridional  rotation  p  are  expressed  in  terms  of  the  displacements  v  and  w 
(see  Figure  A.  1)  as  follows: 


=?r(i+w) 


(A.  2) 


60  =  ^-  (v  cot  0  +  w) 


(A.  3) 


K 


K.=2ii 

6  r2 


~  _L_  iL  [JL  f 

"  ri  d0  \  r1  \d0  "  VJ_ 

r  i  f  dw  —A 

pi  id0 "  VJ_ 

rx  [dp  VJ 


(A.  4) 


(A.  5) 


(A.  6) 


where  r^  and  r^  are  the  radii  of  curvature  of  the  shell  in  the  meridional 

and  hoop  directions,  respectively.  Hookes  law  for  an  orthotropic  shell 

with  the  principal  directions  in  the  hoop  and  meridional  directions  takes 
q 

the  form: 


A.  1 


(A.  7) 


C 1 1  C12  0  0 


C12  C22  0  0 


0  0  C33  C34 


.°  0  C34  C44 


The  shell  configurations  to  be  considered  for  the  bulkheads  and  the  tank 
walls  are  listed  below  with  their  corresponding  geometric  parameters 
defined: 

1.  Conical  Shell  (Figure  A.  2) 
a)  General  Case 


0  =  0, 


r  ^  =  oo 


r2  =  0^)  "  S  cot  (f>r 


(A. 8) 


b)  Cylinder  (radius  =  R) 


,  _  TT 

0  -  T 


=  oo 


(A.  9) 


r2  ~  R 


2.  Ellipsoid  (Figure  A.  3) 

a)  General  Case  (a  =  semimajor  axis,  b  =  semiminor  axis) 


‘1 


'2  .  2  A  .  ,  2  2  V 

a  sm  0  +  b  cos  0  I 


(A.  10) 


[2.2  ,2  Z  V 

la  sm  0  +  b  cos  0j 


A.  Z 


b)  Hemisphere  (radius  =  R) 


(A.  11) 


For  the  present  analysis,  the  longitudinal  displacement  u  and  the  radial 
displacement  v  will  be  more  convenient.  They  are  related  to  the  dis¬ 
placements  V  and  w  by  the  transformation 


v  =  -u  sin  0  +  v  cos  0 
w  =  u  cos  0  -1-  v  sin  0 


(A.  12) 


Substitution  of  this  transformation  into  the  strains,  curvature  and  rotation 
yields 


1  f  du  .  dv 

?7  ta?“n*  +  a? 


COS  0 


(A.  13) 


V 

€e  “  r 

(A.  14) 

•ois 

I  H 

-H  1  u 

II 

■o- 

1  /'du  dv  .  \ 

77  (d0  C°S  *  +  d0  Sm  t) 

(A.  15) 

Jr  COS  0 

Ke  =  ; 

1  fdu  ,  *  d.v  •  A 

77  fe cos  * +  d7 sm  *J 

(A.  16) 

l 

p  “  2 
r  _ 

/■  du  dv  .  \ 

(d0  C°S  ^  +  d0  Sm^j 

(A.  17) 

rl 
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Figure  A.  3.  Ellipsoidal  Shell 


APPENDIX  B 


INITIAL  STRESSES  IN  SHELL  ELEMENTS 


In  the  formulation  of  the  shell  stiffness  matrix,  the  effect  of  the  initial 
meridional  stresses  is  included.  Expressions  for  these  stresses  are 
developed  below  for  each  of  the  shell  elements  to  be  considered.  These 
derivations  are  based  on  membrane  theory  which  is  a  good  first  order 
approximation  except  in  very  localized  areas  where  bending  predominates. 


Conical  Element 

The  derivation  of  the  initial  stresses  in  a  conical  element  is  divided 
into  three  steps.  The  notation  is  listed  below  and  shown  in  Figure  B.  1* 


1. 


R- 


ro 

L 


w. 

x 

h. 

1 


w 

e 

h 

e 

Pi 

pe 

±w 


Stress 


lower  radius  of  the  conic 

meridional  angle 

length  of  element 

effective  density  of  interior  fluid 

location  of  interior  fluid  surface 

effective  density  of  exterior  fluid 

location  of  exterior  fluid  surface 

uniform  internal  pressure 

uniform  external  pressure 


reactive  force  at  top  of  conic 


f+  produces  tensile  stress^ 
-  produces  compressive  t 
<  stress  ^ 


due  to  uniform  pressures  and  reactive  force 


N 


§l  Ztyt  sin  c|> 


W  + 


(Pi- 


Pe) 


- 


X 


(B.  1) 
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2.  Stress  due  to  interior  fluid 


Case  1:  h.  L 
x  ‘ 


N,  =  - .  j  *  pr  dr 

4*2  r  sin  4>0  J  ^ 


-l  fh‘ 

s — t—  I  1  w.  (h.  -  x)  (R0  -  x  cot  cb  )  (-d  x  cot  cb  ) 

;in  cb  I  l  i  2  To'  To' 

To  Jx 


w.  cot  4>  r-  1 

- — r — t“  h.  (h.  -  x)  R.  -  5-  (h.  -  x  (R9  +  h.  cot  6  ) 

r  sin  <po  i  i  '2  Z  I  i  y  '  2  i  To' 


+  T  (hi  "  x3)  cot 


x*ch.  L»  (B.  2) 

^  l 


Case  2:  h.  >►  L 

i 


^2  r  sin  (j) 


pr  dr 


o 


w.  cot  <j>  r  1  ?  ? 

=  +  -1.  ■■  -  x  °  h.  (L.  -  x)  R,  -  1  (L  -  X  )  (R,  +  h.  cot  4  ) 

r  sin  cb  1  1  2  2  v  7  '  2  i  To 

o  L 


+  (LT^  -  x^)  cot  <J> 


0<x<L  (B.  3) 


3.  Stress  due  to  exterior  fluid 


Case  1:  h  ■*=  L 

e  ~ 


-1  rhe  , 

N,  =  - : — 7—  I  pr  dr 

+3  r  sln  4  4 


w  cot  cj>  r  i  /  ?  ?\ 

-  e  r—r  -  h  (h  -  x)  R_  -  i  (h  -  X  (R,  +  h  cot  <j>  ) 
r  sm  cpo  ^  e  '  e  7  2  2\e  2  e  To 


"  2  (h 


,  i  /,  3  3  A  .  , 

+  j  -  x  cot  <j>( 


0  <  x<=  (B.  4) 
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Ellipsoidal  Element 


The  expressions  for  the  initial  stress  for  an  ellipsoidal  element  are 
developed  separately  for  the  upright  and  inverted  bulkheads.  The  notation 
for  both  cases  is  listed  below  and  shown  in  Figures  B.  2  and  B.  3. 


a  radius  of  base 


b  height  of  element 


0o  slope  of  meridian 

(Note:  The  semimajor  and  semiminor  axes,  a  and  b, 
can  be  computed  from  a,  b,  and  0  .  ) 


w^  effective  density  of  interior  fluid 

h.  location  of  surface  of  interior  fluid 

i 

w^  effective  density  of  exterior  fluid 

h  location  of  surface  of  exterior  fluid 

e 

Pj.  uniform  internal  pressure 

p^  uniform  external  pressure 


B.  3 


Upright  Bulkhead  (Figure  B.  2) 

1.  Stress  for  exterior  fluid 


Case  1:  h  fib 

e 


where 


Per2  1  fr 

I.  = - ^ —  +  5 - * — x  2nrp  dr 

<j>^  2  2irr  sin  cp  J—  x 


U 


h  +  b  -  b 


=  a  \  /I  - 


p  =  -  wg  (h0  -  x) 


(B.  7) 


x  +  (b  -  b)  =  b  ~\  1  -  V 


This  finally  leads  to 


Pr_  w  H  2  0  2-  /  ^ 

N.  =  -  ZSJ. - 2-x  -f  (r2  -  r2)-  H  3  +  2jl  1  -  H. 

(j)  2  r  sin  <{>  2  ,,  2  e  3  \  <■ 


Per2 
N  = - e  •  — 

<t>,  ^ 


for  r  s  r  c  a 


for  0  r  r 


(B.8) 


Case  2:  h  =*b 

e 


Per2 


2  27rr 


*_  f 

sin  *  JD 


2trrp  dr 


r_>  w  Hr  2  /  2 

e  2 _ e _  e  a  b  (1  r 

2  r  sin  <j>  2  3  i  ^2 


2\3^Z  2 


for  Of  r<a  (B.  9) 


B .  4 


2.  Stress  for  interior  fluid 


N 


Pir2 


1 


2  2 nr  sin 


j: 

y 


27irp  dr 


where 


H. 

=  a  \/l - iy 

b 


or,  after  integration, 

p.r  w. 
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3b  a  / 
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N 


^2  ^ 


for  r  <  r^a 
for  0  <  r  c  r 


Inverted  Bulkhead  (Figure  B.  3) 

1.  Stress  for  interior  fluid 


Pir2 

N .  =  -i-i  + - 

Yj  2  2irr 

for 


or 


where 


•_  r 

sin  $  J0 


27rrp  dr 


r  ^  r  = 


=  a 
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h.  >0 
1 


h.  <0 
i 


x  ^h. 
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(B.  11) 


(B.  12) 
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This  leads  to  the  following  expressions  for  the  stress: 


„  .  *V2  .  wi  K  2  . 

2  r  sin  4>  L  Z 


p.r_  w.  H.  7  2, 

T  _  i  2 _ i  i  -=2  ,  a  b 

2  r  sin  <j>  l_2  3 


2  +  ^  " 7)  7  0=r=F 


2.  Stress  for  lower  fluid 


Per2 


2  2trr 


‘ _ f 

sin  <j>  J 


2Trrp  dr 


(B.  13) 


r  ■«  r  =  a 


n 


h  +  b 
e 


(B.  14) 


where 


p  =  -wg  (x  -  hg) 


or,  after  integration 


p  r0  w  H  9  2  f  2 

N  =  -  Pe  L  +  e _  _e  r2  4-  ^  [l  - 

<j>?  2  r  sin  cj>  L2  3  l  I  3  ^ 


xih 


ijfor  o 


5  r  r 


-p  r  w  H  2  A  =  2\^^  2 

tvT  *e  2  ,  e  e  -2  ,  a  b  h  r  1  a  b 
N4»7  =  +  715ET?  —  r  +—  \l--z)  -  — 


(B.  15) 
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APPENDIX  C 


FLUID  MASS  MATRIX  EXPRESSIONS 


Detailed  expressions  are  given  below  for  evaluating  Equation  (4.4)  to 
obtain  the  fluid  component  mass  matrix.  Three  cases  are  involved  de¬ 
pending  upon  whether  the  upper  and  lower  tank  bulkheads  are  convex  up¬ 
ward  or  convex  downward  (see  Section  2.0).  A  general  description  of  the 
program  operations  utilized  to  construct  the  mass  matrix  for  the  fluid 
tanks  is  described  in  Appendix  E.  A  brief  flow  chart  followed  by  the 
relevant  equations  in  a  form  suitable  for  a  weighting  matrix  integration 
scheme  is  included  to  provide  the  user  with  a  more  basic  understanding 
of  the  digital  program. 

C.  1  CASE  I.  UPPER  BULKHEAD  CONVEX  UPWARD, 

LOWER  BULKHEAD  CONVEX  DOWNWARD 


Equation  (4.  6)  for  the  general  tank  shell  can  be  expressed  as  the  sum 

of  contributions  to  M,  from  fluid  in  three  sections  of  the  tank,  as  illus - 

b 

trated  in  Figure  C.l  for  Case  I.  One  of  these  sections  is  that  located  above 
the  base  of  the  upper  bulkhead  (shell  1).  A  second  lies  below  the  top  of 
lower  bulkhead  (shell  2).  The  third  lies  between  the  two  bulkheads.  In 
general,  the  motion  of  the  fluid  in  any  section  may  be  affected  by  the 
generalized  coordinate  distortions  associated  with  all  three  tank  shell 
components.  Thus,  Equation  (4.6)  is  expressed  as 


{r(m)°<m) 


(C.l) 


where  m^  =1,  2  or  3,  depending  upon  whether  the  fluid  surface  lies 
within  the  range  of  tank  section  (1),  (2)  or  (3),  as  illustrated  in  Figure  C.L 
H(m)  is  the  depth  of  fluid  within  tank  section  (m),  is  the  radius  of 

shell  (m)  and,  similarly,  u(x)  and  r)  represent  the  fluid  motions 


within  tank  section  m.  The  specific  form  of  the  column  matrix  /u,  %(xfV 

V.  (m)  v 

m  terms  of  the  component  shell  generalized  coordinates  is 


C.  2 


and 


0 


are  the  generalized  coordinate  displacement  vectors  for  tank  shell  com¬ 
ponent  m,  and  is  the  meridional  angle  for  tank  shell  component  m, 

and  U  and  V  are  the  number  of  longitudinal  and  radial  generalized  coordi¬ 
nates,  respectively,  in  shell  component  m. 

The  column  matrix  r)}  i*1  terms  of  the  component  shell 

generalized  coordinates  is 


C.  3 


0 


r(2)^(2)(x’r) 


2  cot  ?(z)fu)\  \ 

(  'l  }  -— - ^'<4fLS  ,(x))  (C.  6) 

l2  (2)  J 


C.  2  CASE  II.  UPPER  BULKHEAD  CONVEX  DOWNWARD, 

LOWER  BULKHEAD  CONVEX  DOWNWARD 

The  Case  II  configuration  is  illustrated  in  Figure  C.  2  which  also 
defines  the  tank  sections  m  for  this  case.  Equation  (C.  1)  for  the  fluid  in 
tank  section  (1)  is  modified  as  follows.  The  fluid  motion  ^u(x)^  in  tank 
section  (1)  is  given  by 


C.  4 


(C.  8) 


The  fluid  motion  ^v(x,  r)^  in  tank  section  (1)  is  defined  by 


C.  5 


Upon  substitution  of  Equations  (C.  1)  and  (C.  2)  into  Equation  (4.  1)  and 
integration  with  respect  to  r  between  the  limits  of  and  r(2)’  one 

obtains  the  following  expression  for  the  m  =  1  portion  of  Equation  (C.  1). 


(C.  10) 


C.  6 


where 


The  expressions  for  the  m  =  2  and  m  =  3  portions  of  Equation  (C.  1) 
remain  unchanged.  In  summary  for  evaluating  the  equivalent  of  Equation 
(C.  1)  for  Case  II  tank  configuration,  use  Equations  (C.  10),  (C.  11),  (C.  12) 
and  (C.  13)  for  m  =  1,  and  Equations  (C.l),  (C.3),  (C.4),  (C  6)  and  (C.  7) 
for  m  =  2  and  m  =  3 . 

C.3  CASE  III.  UPPER  BULKHEAD  CONVEX  UPWARD,  LOWER 
BULKHEAD  CONVEX  UPWARD 


The  Case  III  configuration  is  illustrated  in  Figure  C.  3  which  also 
defines  the  tank  sections  m  for  this  case.  Equation  (C.  1)  for  the  fluid 
in  tank  section  (3)  is  modified  as  follows.  The  fluid  motion  ^u(x)^  in  tank 
section  (3)  is  given  by 

r  a  1 


C.  8 


Upon  substitution  of  Equations  (C.  14)  and  (C.  15)  into  Equation  (4.  1)  and 
integration  with  respect  to  £  between  the  limits  of  r^  and  r^,  one 
obtains  the  following  expression  for  the  m  =  3  portion  of  Equation  (C,  1). 


C.  10 


C.  11 


Equation  (C.  1)  remains  valid  for  m  =  1  and  m  =  2,  however.  Equations  (C.2) 
and(C.3),  respectively,  become 


C.  12 


0 


In  summary  for  evaluating  the  equivalent  of  Equation  (C.  1)  for  Case  III 
tank  configuration,  use  Equations  (C.  1),  (C.  20)  and  (C.  5)  for  m  =  1, 
use  Equations  (C.  1)  with  limits  of  integration 


Equations  (C.  21)  and  (C.  6)  for  m  =  2,  and  use  Equations  (C.  16),  (C.  14) 
and  (C.  15)  for  m  =  3. 


C.  13 


C.  15 


APPENDIX  D 


CALCULATION  OF  THE  SHELL  STIFFNESS  AND  MASS  MATRICES 

In  this  section  the  equations  for  the  shell  stiffness  and  mass  matrices 
presented  in  Section  3.0  are  rewritten  in  a  form  more  suitable  for  digital 
programming.  As  previously  mentioned,  integration  of  these  equations  is 
accomplished  using  a  sixteen  (16)  point  Gaussian  weighting  matrix  integra¬ 
tion  schemed 

These  operations  may  be  performed  in  a  straightforward  fashion  using 
matrix  operations.  The  equations  in  this  section  are,  therefore,  written 
in  matrix  notation  and  are  identical  to  the  equations  which  are  coded  into 
the  digital  program.  For  a  more  complete  understanding  of  the  steps  re¬ 
quired  in  the  construction  of  the  shell  stiffness  and  mass  matrices,  a 
functional  flow  chart  depicting  the  program  operations  is  presented  in 
Section  D.  1.  The  weighting  matrix  coefficients  are  given  in  Section  D.  Z 
and  other  detailed  expressions  are  provided  in  Sections  D.3,  D.  4  and  D.  5 

for  elements  KM,  M  and  for  the  initial  stresses  N  . 

0 


D.  1 


D.  1  FUNCTIONAL  FLOW  CHART  FOR  CALCULATION  OF  SHELL 
STIFFNESS  AND  MASS  MATRICES 


1  RESET  INDEX 


ASSEMBLE  SYSTEM 
COORDINATE 
VECTOR  AND  STORE 
f  COORD.  "\ 

{C)a  INPUT  / 
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G  SEE 

10  Section 
„  2  D.Z  FOR 
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CALL 

BASIC  DATA  FOR 
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TEST  SIGN 
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M  =  1,  2,  13 

NOTE:  SEE  SECTION  D.  3 


CONICAL 

SHELL 


ELLIPSOIDAL  SHELL 


ASSEMBLE 


[*w]  (0  +  V)x(0  +  V)’  +  [Kw] 


“GO  TO  NEXT  PAGE 
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D.  1  FUNCTIONAL  FLOW  CHART  FOR  CALCULATION  OF  SHELL 
STIFFNESS  AND  MASS  MATRICES  (Continued) 
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NOTE:  SEE  SECTION  D.  3 


RESET  INDEX 
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D.  2  GAUSSIAN  WEIGHTING  MATRIX  TABLES 


DATA  FOR 

[*Si4] 

DATA 

wgJ 

G1 

0.00529 

95325 

W1 

0. 01357 

62297 

G2 

0. 02771 

24885 

W2 

0.03112 

67620 

G3 

0.06718 

43988 

W3 

0.04757 

92558 

G4 

0.  12229 
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W4 

0. 06231 

44857 

G5 
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W5 

0.07479 

79944 

G6 

0. 27099 
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W6 

0.08457 

82597 

G7 

0.35919 
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W7 

0.09130 

17075 

G8 

0.45249 

37451 

W8 
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53052 

G9 
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W9 
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53052 

G10 
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W14 
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G15 

0. 97228 
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W15 

0.03112 

67620 

G16 

0.99470 

04675 

W16 

0.01357 

62297 

D.  3  ASSEMBLY'  OF  [rKMF]  FOR  M  =  1 ,  2,  •••,  13: 

The  equations  for  KM  and  F  are  presented  below  and  are  evaluated 
at  the  points  £  =  G^,  G^,  *  ■  •  ,  G^. 
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sin^  0  +  "y  |  sin  0  +  cos  0j 
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N  cos 
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KlO  -  K10  -  -C|  j  sin  0  cos  0  -  C33  — — ^  |cos  0  -  -pi-  sin  0^  |sin  0  +  cos  ^ 

_  sin  0  cos  0  f  ,  .  rl  A 

"  C34 - r.'r1  ^in0+  — cos  0j 


'34  r 
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^  sin  0  cos  0  ,  .  . 

t  +  N  sin  0  cos  0 

r  ^ 


in  =  K1 1  =  -C  -jr-  ^sin  0  +  cos  0^  sin  0  +  -~^n  ^cos — $_ 


KT2  =  K12 


=  C33(cos  0  -  -p  sin 
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sin  o  cos 


"34  r 


K13  =  K13  =  CL-  sin  0  cos  0 


r  Z 

1  -3(k  -  1)  sin  0  cos  0 

■  w  - 1)  iz?  ♦ ;  ii 


M  =  2tt v  rtF 

1  a 


F  =  _ — —  £or  conic 
sin  o 

o 


=  r^0Q  for  convex  upward  ellipsoidal  bulkhead 
=  r^(0Q  -  tt)  for  convex  downward  ellipsoidal  bulkhead 


D.  4  INITIAL  STRESS  FOR  CONIC 
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D.  5  INITIAL  STRESS  FOR  ELLIPSOIDAL  BULKHEAD 
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APPENDIX  E 


CALCULATION  OF  FLUID  MASS  MATRIX 

The  construction  of  the  fluid  mass  matrix  presented  in  Section  4.  0 
for  three  tank  configurations  differs  from  the  construction  of  the  shell 
inertial  and  stiffness  characteristics  insofar  as  the  fluid  equations  involve 
double  integrations.  As  mentioned  in  Section  4.0,  these  integrations  are 
performed  using  a  Lagrangian  integration  scheme.  Although  this  tech¬ 
nique  is  similar  to  the  Gaussian  method  of  integration  used  in  the  compu¬ 
tation  of  the  shell  matrices  described  in  Appendix  D,  only  the  Lagrangian 
technique  is  adaptable  to  double  integration.  The  digital  program  utilizes 
a  21 -point  Lagrangian  weighting  matrix  which  is  constructed  by  combining 
two  11 -point  matrices  tabulated  in  the  above  reference.  The  increased 
number  of  points  provides  additional  numerical  accuracy.  In  the  following 
pages,  the  fluid  equations  are  rewritten  to  accommodate  the  integrations 
using  the  weighting  matrix  technique.  A  functional  flow  chart  outlining 
the  complete  mass  matrix  calculation  is  provided. 


The  notation  of  this  appendix  corresponds  to  the  notation  utilized  in 
the  digital  program. 
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See  Functional  Relationship  £a2(£al)  defined  below. 
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See  functional  relationship  £^2  (£  )  defined  below 
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Functional  Relationship  £a2(£al) 
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Functional  Relationship 


ra2  =  |ry^b-b 


1  -n  1  - 


a2 


where 


r  =  f<Sa3> 


4  .  2  ^ 

a  sm  0 

“2  ~ 2  71  T 


(a  sin  0  +  b  cos  0) 


and 


£  =-£- 

^a3  ^ 


E.  8 


E.  3  CASE  I  /al  -  Convex  Up  Ellipsoid;  a2  -  ConicT 
a3  JIA:  Convex  Down  Ellipsoid^ 


iIB:  Conic 


E,  3.  1  Integral  Equation  Relationships 


re  ii 

^10 


d^al 


E.9 


E.  11 


E.  12 


r 


E.  3.  Z  Integration  Limit  Data 
Fluid  Surface  in  Section  1 


*n=K)al 

^20  =  ^30  =  ° 

^10  =  ^21  =  ^31  = 
Vi  =  ^2  =  73  =  *Y 


Note: 


<p  =  sin 


1 


2.2  -2.,  2 
a  k  -  r  (k  -  1) 


T7T 


k  =  a/b 


if  (Hj  +  b  -  E)  =»  0  ,  0  x  0  «c  90°,  if  <=  0,  90°  «=  <j>  <=  180° 


Fluid  Surface  in  Section  2 


$10  “  $11  ~  ^31  1 

^20  =  $30  =  ° 

V° 

yz  =  y3  =  y 

Fluid  Surface  in  Section  3 


$10  =  $11  =  1 

$20  =  $21  =  $30  =  ° 

^i  =  72  =  o 

y3  =  y 


E.  14 


Note: 


0  =  sin 


k  =  a/  b 


-[a2k2-?V-  I)]' 


{  .  (b  -  H3)2 


E.3.3  Specific  Matrix  Data 


Note:  Quantities  <(  ) 


are  evaluated  at 


am,  1 


points  £am  j  (i  =  1,2,3,...  21)  where 

t  ,  =  £  n;  £  7  =  (0.  95  £  n  +  0.  05  ?  .); 

bam,  1  bm0  bam,  2  bm0  bml 

?  .  =  (0.90?  n  +  0.  10  £  .);  ?  .  =  (0.  85  ?  n+0.  15?  _); 

bam,  3  v  bm0  bml  bam,  4  v  bm0  bml' 

7n  =  (0.  05  £  A  +  0.  95  £  J;  £  71  =  g  . 

bam,  20  bm0  bml  bam,  21  bml 


(^ll  "  ^10) 


K2.  = 


i  +4it^2  ^L^a2  \  (  2/  /  '^21 

'*2  I  ,  ; 


IAJK3i=-4^3(,o-.W{li^ 


£?r>) 


^31  "  ^30^ 


IB:  <  K3 


5.  -  +4tt^3  (|L|)a3  <  f  21  >  '  ^3  1 

l^r  'a3  J  ,  . 


E.  15 


E.  16 


(kU2)  =  -(|L|cot0o^|(r) 


a2  J  .  '  (^21  “  ^20)  ; 

Ja.Z,  1 


KUZ  =  KUZ  =  KV2  =  KV2  =  0 


4  •  (S21  -  e20) 

[Kvi\  ~~  +  (*°  -  ”)a3(lr  •  rl  cos  )a3  .  •  <«31  -  «30> 


KV3j.  =  +  ^o  -  ,rJa3\|r  ■  rl  sin  »>a3/a3_  .  ’  '631  -  ho] 


(kU3)  =  -(|L|cot0o)  |(r) 

'  'i  a3 1, 


a3  [  _  .  ’  ^31  "  ^30^ 
)  a3,  1 


(KV3)i  =  .,iLDa3(ir)a3]  •  ,e31  -  e30) 


H^Wa.,3 

KhN^Ri 


KV5 j.  =  +  [<r>a2 Ja2>  .  '  KU5  =  KU5  =  KV5  =  KV5  =  0 


E.  17 


< 


E.  18 


E.  4  CASE  II  /a  1  -  Convex  Down  Ellipsoid;  a2  -  Conic;\ 
fllA:  Convex  Down  Ellipsoid*! 

\  \IIB:  Conic  J  j 


E.  4.  1  Integral  Equation  Relationships 


m- 

V  31  *  10  va2  '  ralj  l 


«,i  ('■  i,  /('u  -  ■!,) 


u(x)) 


\T\l  -  ral)  A  A 

- 2 -  u(x7  d^a  1 


1  T 


+4*V2  (|L|)a2  V7-“(xy 


a2 


>|c 

Identical  with  Case  I. 


E.  19 


E.  21 


E.4.2  Integration  Limit  Data 
Fluid  Surface  in  Section  1 


E.  4.  3  Specific  Matrix  Data 


Same  as  for  Case  I 


E.  23 


(kv4)  -  (1) 

x  'l 


‘Same  as  for  Case  I 


E.  24 


Same  as  for  Case  I 


E.  25 


(K5)i  =  -Z  j(i 


(K5).  =  0 


nA:{(K6)t- 


IIB:  {(K6).  =  -2(cot0o) 


a3,  i 


Same  as  for  Case 


I 


E.  26 


E.5  CASE  III  (al  -  Convex  Up  Ellipsoid;  a2  -  Conic; 
a3  -  Convex  Up  Ellipsoid) 


E.  5.  1  Integral  Equation  Relationships 


a3 


E.  27 


E.  29 


0 


E.  30 


E.  5.  2  Integration  Limit  Data 
Fluid  Surface  in  Section  1 


^irK, 

*20=«b»a3/<|L|>a3 

«31'° 

£  10  =  ^21  =  ^  30  =  1 

Vi  =  V2  =  y3  =  v 

Fluid  Surface  in  Section  2 

*10  =  *11  =  ^  30  =  1 
e2i  =  (H/|LDa2 
g2o  =  (lbl)a3/(IM)a2 

e3i  =  ° 

Vl  =  ° 

v2  =  v3  =  \ 

Fluid  Surface  in  Section  3 

£io =  =  1 

&20  =  =(53/*I-l)a2 

^30  =  1 

«31-(W»o)l3 


E.  32 


Yl  =  Y2  =  0 

y3  =  y 

E.  5.  3  Specific  Matrix  Data 


(C2).  =  (C3).  =  0* 


'n'Y? 

(C4).  =  +-2-(lLl)a2\1}a2  i  '  (£21  -  iz0) 

’  '  ^31  "  ^ 30^ 

a3,  i 

(C6).  =  - WV3 (^0)a3 p2< ra3 ’  ra2»  ’  (rl  sin  *^3)3  .  '  '«31  '  «30> 

(C7).  =  +"V3(^0)a 3{F1(ra2-  ra3*  '  (rl  sin  ♦)  }  '  <«31  '  «30> 

^a3,i 


(C5).  =  -^3^0)a3^F1(ra3,  ra2)  •  sin  ^ 


❖ 

Same  as  for  Case  I 


E.  33 


<KUi>i  =  +KJ  A  lr  *  ri cos  , )  '  (^ii  -  ^io> 

aU  al4i,i 


(KVl)i  =  +  J  Mr.rlSin0  l  *  <^11  "  ^10> 
al(-  al,  i 


(KU2).  =  -(|L|cot0o)a2|(r)a2Ja2  _  •  (£21  -  t2Q) 


<KV2)i=-<lLl>a2(<r)a2Li*^21-^0) 


* 


(KU2).  =  (KV2)  =  0* 


<KU2>i  =+(cot«,0)a2(«.0)aj|(r)a2(r1  sin  .  ■  (g31  -  l30> 


*  ^31  "  ^30^ 


<^2>i  =  +(^o)a3{(r)a2(rl  Sin  *) 

<KU3>i  =  'Wa3{(r  •  rl  COS  *)a 

(KV3).  =  •  rj  sin  ^ J  *  d3l  '  t3Q) 

3.3)  1 


a3,  i 


a3,  i 


^31  "  ^30* 


(KU4).  =  +<  (r  cot  *)al 


al ,  i 


(KV4>i  =  +  (r>al| 


al ,  i 


(KU5)^  =  (KV5)  =  0" 


❖ 


Same  as  for  Case  I 


E.  34 


► 


* 


r 


(KU5).  =+(cot0o)  <  ( r) 


a2l 


a2 


a2,  i 


(KV5).  =  + 


(KU^)i  =  + 


j(cot  <p) 


a3,  i 


(KV6).  =  +1 


(KU5).  =  +  j(cot  0o)  } 

\  a3y a3,  i 


2  cot  <p 
r 


al ,  i 


(KV5).  =  +1 


(K4)i  = 


(K5).  =  0 


(K5).  =-2(cot0  ) 

1  °  a2Ur7a2 


a2,  i 


(k5).  = 


<K6>i  = 


-2(cot0o) 


1 


a2  1  (ra2  -  ra3 ) 


a3,  i 


“2  ((cot  »)a3  •  J  2  1  7-, 

L  (ra2  “  ra3 j 


a3,  i 


W  ’ 

Same  as  for  Case  I 


E.  35 


E.  6  LAGRANGIAN  WEIGHTING  MATRIX  TABLES 


DATA  FOR 

M 

DATA*  FOR  [w] 

21x21 

wi 

0. 01341 

70742 

W1 

0. 01341 

70742 

W2 

0. 08876 

79707 

W2 

0. 08876 

79707 

w3 

.  -0.04052 

17853 

w, 

-> 

-0.04052 

17853 

w4 

0. 22747 

31442 

W4 

0. 22747 

31442 

w5 

-0.  21757 

75613 

W5 

-0. 21757 

75613 

W6 

0. 35688 

23152 

W6 

0. 35688 

23152 

W? 

-0. 21757 

75613 

w7 

-0.21757 

75613 

W 

8 

0. 22747 

31442 

W8 

0. 22747 

31442 

W 

9 

-0. 04052 

17853 

w9 

-0. 04052 

17853 

W10 

0. 08876 

79707 

W10 

0. 08876 

79707 

W11 

0. 02683 

41484 

W 

11 

0. 02683 

41484 

W12 

0. 08876 

79707 

w 

w  12 

0. 08876 

79707 

W13 

-0. 04052 

17853 

W 

W  13 

-0.04052 

17853 

W14 

0. 22747 

31442 

W 

W  14 

0. 22747 

31442 

W15 

-0.  21757 

75613 

W15 

-0. 21757 

75613 

W16 

0. 35688 

23152 

W1 6 

0.  35688 

23152 

W17 

-0. 21757 

75613 

W17 

-0.21757 

75613 

W18 

0. 22747 

31442 

W18 

0.  22747 

31442 

W19 

-0. 04052 

17853 

W19 

-0.  04052 

17853 

W20 

0.  08876 

79707 

W20 

0. 08876 

79707 

W21 

0. 01341 

70742 

W21 

0. 01341 

70742 

This  column  is  identical  for  all  21  columns  of  the  matrix. 
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.06849  0.  08623  0.08528  0.08554  0.08541  0.08551  0.08540  0.08560  0.08498  0.08876  0.08876  0.08876  0.08876  0.08876  0.08876  0.08876  0.08876  0.08876  0.08876  0.08876 

51420  68393  88596  06980  85560  29870  09862  69692  04181  79707  79707  79707  79707  79707  79707  79707  79707  79707  79707  79707 
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.01928  -0.01603  -0.01713  -0.01651  -0.01706  -0.01629  -0.01822  0.00256  0.05239  -0. 04052 t 

76386  82195  27364  69152  38131  05844  10802  40532  81078  17853; 


APPENDIX  F 


SAMPLE  INPUT  DATA 

Sample  input  data  for  a  typical  one  stage  launch  vehicle  (Figure  F.  1) 
is  presented  to  illustrate  the  basic  data  requirements.  For  simplicity 
all  structural  components  are  constructed  of  aluminum.  The  vehicle 
oxidizer  and  fuel  tanks  are  both  assumed  to  be  simple  tanks  according 
to  the  terminology  used  in  Section  2. 0.  To  demonstrate  the  capability  of 
handling  stringers  and  ribs,  two  shell  sections  are  given  orthotropic 
properties  as  indicated  in  Figure  F.  1. 

The  physical  model  is  first  subdivided  into  a  consistent  set  of  shell, 
fluid  and  mass-spring  components  as  shown  in  Figure  F.  2.  In  this  example 
the  vehicle  is  represented  by  eleven  (11)  shell  components,  two  (2)  fluid 
components  and  four  (4)  spring-mass  components  to  account  for  the  pay- 
load,  engine  and  equipment.  The  displacement  coordinate  locations  are 
then  selected  and  numbered  according  to  the  requirements  discussed  in 
Section  8.0.  For  this  sample  problem  the  vehicle  is  assumed  to  be  unsup¬ 
ported,  i.  e.,  no  fixed  coordinates  are  specified. 

The  input  data  load  sheets  have  been  prepared  for  the  launch  vehicle 
illustrated  in  Figure  F.  1  and  are  included  in  this  appendix  to  illustrate 
the  input  format.  The  system  data  appears  first  followed  in  order  by 
the  data  specifications  for  each  of  the  shell,  fluid  and  spring-mass 
components. 


F.l 


1400  POUNDS 


Figure  F.  1.  Typical  One-Stage  Launch  Vehicle 
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1 


LEGEND 


Q  SHELL  COMPONENTS 

A  FLUID  COMPONENTS 
□spring-mass  COMPONENTS 


Figure  F.  2.  System  Components  and  Coordinates 
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